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DATA REDUCTION PROGRAMS 
FOR A LASER RADAR SYSTEM 
By 

F. F. Dadavi 1 and G.E. Copeland 2 


SUMMARY 

The purpose of this report is to sunmarize the state of work in the de- 
velopment of software routines, that are used in the pulsed C0 2 laser radar 
system currently under development at NASA-LARC. 

This work was performed under research grant NAG1-1, during the period 
of August 1 to October 15, 1983. 


INTRODUCTION 

This report gives the listing and description of software routines 
which were used to analyze the analog data obtained from LIDAR - system. 

All routines are written in FORTRAN - IV on a HP - 1000/F mini -computer 
which serves as the heart of the data acquisition system for the LIDAR 
program. This particular system has 128 kilobytes of highspeed memory and 
is equipped with a Vector Instruction Set (VIS) firmware package, which is 
used in all the routines, to handle quick execution of different long loops. 
The system handles floating point arithmetic in hardware in order to enhance 
the speed of execution. This computer is a 2177 C/F series version of HP - 
1000 RTE-IVB data acquisition computer system which is designed for real 
time data capture/analysis and disk/tape mass storage environment. 


Research Associate, Department of Physics, School of Sciences and Health 
Professions, Old Dominion University, Norfolk, Virginia. 

2 Associate Professor, Department of Physics, School of Sciences and Health 
Professions, Old Dominion University, Norfolk, Virginia. 


The main prograns in this report are LIDAR, TEMPE, and DPIOT. Each 
program is composed of a main routine which either takes care of Initializ- 
ing all the default values, or does the proper calling sequence to other 
segmented routines to handle data acquisition, storage, retrieval, computa- 
tion, smoothing and plotting. 

The need for a segmented approach to programming Is due to memory limi- 
tation of the system in use. The benefit of a segmented program is that it 
does not have to reside in the main memory until it is needed. Upon calling 
a swapping between disk and main memory a segmented routine can do its func- 
tion and be replaced by another segment. This modular approach to program- 
ming eliminates the need for large memory requiranent. 

The remainder of this report will give a brief description of each pro- 
gran, with commented listing of the routine Itself at the end. Included are 
some of the results obtained. 


PROGRAM LIDAR 

LIDAR, which was partly developed prior to the start of this work, is 
the main routine for data handling. It accepts digitized data from an A/D 
converter, (which is the Biomation 1010 digital waveform recorder) . This 
converter provides 10 bit resolution and sends data to the computer in a bit 
parallel/word serial fashion. During this work, the program and its seg- 
ments were modified to make them capable of storing data on a 9-track mag- 
netic tape recorder. 

The main program LIDAR is written in an interactive fashion, so that 
the operator can choose the preset default values or enter new values 
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through the system console. The operator Is given the option of plotting 
the results on a graphics terr!nal or hardcopy plotter in order to preview 
the data and decide upon further options such as storage/ dimp, nimber of 
curves/number of frames. Figure 1 is a typical output of this program and 
Its segments for a particular run and it shows the attenuation of returned 
signal versus amplitude. Four curves were plotted within one frame and the 
time of the day is written on the upper right corner of graph for each 
curve. Included is the date and nimber of averaged shots to obtain the da- 
ta. 

The segmented routines for LIOAR are BISEG, TRASH, POSTP and PTLID. 

The function of these segments briefly are: 

BISEG reads data from Biomation and averages NAVG shots, each 
consisting of the first 512 points of Biomation memory. The 
second 512 words are used to determine a baseline level. 

TRASH, which is written for the sole purpose of writing a 
trash file at the end of a good series of data on the magnetic tape, 
performs all of its functions through a series of EXEC calls. A 
trash file is usually needed so that the operator would know when 
the last set of data was taken on the previous run. 

POSTP handles all the calls to magnetic tape recorder and gen- 
erates mechanical motion upon receiving instruction from operator, 
they include forwarding, backward, skipping of data stored on the 
tape. This is done by reading the heading of each data file and 
printing it on the system console, so that the operator can decide 
where the end of file is. 
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PTL ID Is the general purpose plotting package for plotting 
the analyzed data on the plotter/graphic terminal. The 
plotting is done through many calls to the HP - plot package. 

PROGRAM TEMPE 

TEMPE is written for the purpose of fast smoothing the stored data 
on the magnetic tape. The short smoothing technique is applied to all the 
data points and the result is plotted versus raw data. This program was 
written when analysis of the stored data was at an early stage and the 
degree of smoothing was unknown. The program takes advantage of three and 
five points running averaging techniques to smooth the data. The main pre- 
program TEMPE has the same modular structure as LIDAR and in many ways is 
similar to it. It is written in the interactive fashion and allows the 
operator to select the default values for various plots or input new 
val ues. 

The segmented routines for TEMPE are FNTEM and PTTEM. The function of 
these segments briefly are: 

FNTEM takes care of operating the tape recorder and positions 
the tape such that the headings of all the files are printed on 
the system console. The operator has the option of skipping fow- 
ward or backward any nunber of files. The data is then read back 
to TEMPE and quick smoothing is performed on it. 

PTTEM is the general purpose plotting package for plotting 
the smoothed data versus raw data for the sake of comparison. 

The routine has the same structure as PTLID in calling HP plot 
package. 
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PROGRAM OPLOT 

OPLOT Is written for the purpose of applying general smoothing tech- 
niques to the stored data on the tape. The three smoothing routines are 
TRIANGULAR, RECTANGULAR and GAUSSIAN. Each one may be applied to the data 
and the result Is plotted along with an analytical fit to the smoothed 
curve. 

The main program DPLOT calls segment REDUC to properly position the 
tape and gives the operator the option of selecting any one of the three 
smoothing routines which are included in the subroutine SMOTH. After 
performing the required task, control Is returned back to the main program 
DPLOT. At this point the program exits to segment PTRED for plotting the 
smoothed curve and, finally, DAFIT calculates a theoretical fit to the data 
and plots it through PTRED. 

Figures 2(a), 2(b) and 2(c) are typical results obtained. A theoret- 
ical fit is plotted along with a smoothed data curve. The vertical lines on 
the upper x-axis shows the fitting region and information concerning smooth- 
ing and fitting is given on the upper right corner. 

The segmented routines for DPLOT are REDUC, DAF12 and PTRFD. The func- 
tions of these segments briefly are: 

REDUC handles the operation on the tape recorder and 
positions the tape such that the headings of all the files are 

printed on the system console. The operator has the option of 
skipping forward or backward any nunber of files. Once this is 
done a smoothing technique is applied to all the data points. 

DAF12 performs a theoretical fit to the data points based 
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BLOCKED B: 5B 



Figure 2(a). A rectangular fit to data. 
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Figure 2(b). A triangular fit to data. 
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Figure 2(c). A gaussian fit to data. 






on the fitting distance which Is Inputted by the operator. 

PTRED Is the general purpose plotting package for plotting 
the smoothed data versus the fit for the sake of comparison. The 
routine has the same structure as PTLID in calling HP plot pack- 
age. 
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PROGRAM LIDflR 

COMMON/XNY/X ( 5 1 2 ) , Y ( 5 1 2 ) , I TOUT, NAVG , I FRAN, I OPEN, KFI LE 
COMMON/GRAPH/I GCBC192) ,1 PRAM ( 5 ) , 1 0 BUF ( 20 ) , I DCB C 1 4 4 ) , I F I LE ( 3 ) 
XMIN, XMAX, YMIN, YMRX, NCRV, ICRV, VXMIN, VXMAX, VYMIN, VYMAX, NPLT 
* , NPLO T , v ALI B , XLO OK , I E NTR Y , MTO PT , MQ N , I DAY , I YE RR , I MG UR , MI NUT , I AM 
COMMON I NSTN1 , 1 NSTN2 , 1 NSTN3 , 1 NSTN4 , 1 NSTN5 
DI MENSI ON VPNOM ( 4 , 4 ) , I DAYS ( 1 2 ) , I TI ME ( 5 ) 

INTEGER BI SEG ( 3 ) , PTLI DO), TRASH C 3 ) , POSTP ( 3 ) 

EQUIVALENCE (IPRRMC 1 ) , LU) , ( I PRAM ( 2 ) , I D ) , < I PRAM ( 3 ) , NAVGDF ) 

DATA VPNOM/1 0. , 70. , 57 . , 9? . , 1 0, , 70 . , 7 . , 4? . , 05 . , 1 45 . , 57 . , 9? . , 

•05. , 145. ,7. ,47.x 

/ 

/ 

/ 

/ 


BATA 

DATA 

DATR 

DATA 


BI SEG/2HBI , 2HSE , 2HG 
PTLI D/2HPT, 2HLI , 2HD 
TRRSH/2HTR, 2 HAS, 2HH 
P0STP/2HP0, 2HST , 2HP 


C 

C 


DATA IDRYS/31 ,28,31,30,31,30,31,31,30,31,30,31/ 


THIS CALL PICKS UP THE 
CALL RMPAR C I PRAM ) 


LU AND ID OF THE PLOTTING DEVICE 


600 


FORMAT C" 
. •• 


OR",/, 
PLOTTER",/, 




501 FQRMRT( M 


THIS PRORGAM IS THE MAIN PROGRAM FOR SETTING UP PLOTS",/ 
IT READS IN THE PLOT DEVICE LU AND ID THRU RMPRR. 

I.E. RU, LIDflR, 1,1 PLOTS ON THE CRT 
RU, LI DAR, 32,2 PLOTS ON THE 
RU, LIDflR, 33, 2 PLOTS ON JIM'S PLOTTER",/) 

QUESTIONS WILL BE AKSED CONCERNING HOW MANY PLOTS",/ 


C 

C 


" SCALING 
WRITE (1,600) 

WRI rt£( 1 , sol ) 


OF EACH (INDIVIDUALLY) AND PEN NUMBER",/) 


THIS ASSURES THAT R VALID LU/ID COMBINATION IS USED 
IFCLU.LE.G .OR. ID.LE.O) GO TO 90 
IF(ID.EQ.I) LU-1 

IF (ID. EQ,2. AND. (LU.NE . 32 . AND. LU. NE. 33 ) ) GO TO 90 
IF(ID. GT. 2) GO TO 90 


300 ‘ r ORMATC 1 X," HOW MHNY FRAMES DO YOU 


400 


WANT?",/," DEFAULT: 


4" ) 


'?« / 

> ' > 


FORMAT!" INPUT VIEWPORT COORD.: VXMIN, VXMAX, VYMIN, VYMAX" , /, 

FOR FRAME", 13,/," DEFAULT VALUES: ", 3 (F3 . 0, ","), F3 . 0) 

401 FORMAT (" INPUT WINDOW COORD. XMIN, XMAX, YMIN, YMAX" , /, 

DEFAULT VALUES: 0,500,-100,900") 

402 FORMAT (" HOW MANY CURVES ON FRAME", 1 2," 

-"NOTE 0 IS AN ABORT",/. 

-" 4 IS THE DEFAULT" ) 

403 FORMAT (//I OC'#" ), /"INPUT SAMPLE INTERVAL IN MICROSEC" 

DEFAULT VALUE: .2" ) 

404 FORMAT (//I OC" #"),/» INPUT 1 FOR HORIZONTAL OR 2 FOR VERTICAL" , /, 

-" DEFAULT VRLUE: 1 (HORIZONTAL)") 




C 

C 


501 FORMAT (" INPUT 
-" NOTE --- 

^ii 


THE 


AVERAGE" 






502 


FORMAT(//" 

u 


503 FORMRT(" 

. ii 


NUMBER OF SHOTS YOU WANT TO 
0 IS A #* CURVE ** ABORT" 

DEFAULT:" ,14) 

****#*#***H*######**#***#*»******i****#******#*#*" 

* PLACE PEN IN HOLDER #»,/, 

#**#*#***********•**#***'*******************•*****<' , //) 




504 


505 


INPUT DEFAULT VALUE FOR NO, 
DEFAULT: 32" ) 

DO YOU WANT TO STORE THESE 
DEFAULT: N" ) 

FQRMAT(R1 ) 


FORMATC** 
. ii 


OF SHOTS TO BE AVERAGED" , /, 
DATR ON MAGNETIC TAPE?",/, 
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0065 

0066 
0067 
0060 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 
0000 
0001 
0082 

0003 

0004 

0005 
0086 
0007 

ooee 

0009 

0090 

0091 

0092 

0093 

0094 

0095 
Ousts 
0097 
0090 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 
0107 
0100 

0109 

01 10 
01 1 1 
0112 
01 .3 
0114' 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 
0127 
0120 

0129 

0130 


C 

C 


C 


C 


C 


«.n 

-ii 


506 FORMRTC //•’ IS THIS A #4* BLOCKED ### CASE?”,/, DEFAULT: N" ) 

507 FORMRTC" PRESS ## RETURN ## TO CONTI NUE” ) 

508 FQRMATC///2X,30C#” ),/, 

" IF THIS IS THE FIRST RUN OF THE DRY FOR WHICH " ,/, 

YOU RRE PLANNING TO STORE DATA ON TAPE , YOU PROBABLY WANT”,/, 

TO REPOSITION THE TAPE TO THE END OF THE PREVIOUSLY RECORDED”/, 
BATAo THIS REQUIRES A LITTLE TIME SO BE PATIENT.”,//, 

DO YOU WANT TO REPOSITION THE TAPE? (Y OR N)”,/, ) 

509 FORMAT (" FOR THIS QUESTION THERE IS NO DEFRULT. ANSWER Y OR N”/) 

510 FORMRTC POSITIONING TAPE....”) 

MTCNT-0 

IYES-1HY 
N0-1HN 
ASSIGN 101 


TO INSTN1 
ASSIGN 102 TO INSTN2 
ASSIGN 103 TO INSTN3 
ASSIGN 104 TO INSTN4 
ASSIGN 105 TO INSTN5 
IENTRY-0 
NPLOT-O 
ITOUT-50 
WRITE Cl, 503) 

RE AD Cl,#) NAVGDF 
IFCNRVGDF .EQ.O) NAVGDF -32 
WRI TEC 1,300) 

RE AD Cl,#) NPLQT 
IF CNPLOT .EQ.O) NPLOT-4 
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DO 92 NPLT-1 , NPLOT 
VXMIN-G. 

VXMAX-O. 

KNQM-MQDCNPLT, 4 ) 

IFCKNOM .EQ. 0) KNOM-4 

WRITE Cl ,400) NPLT, CVPNOMC JNOM,KNOM) , JNOM-1 ,4 ) 
RERD C 1 , # ) VXMIN, VXMRX, VYMI N, VYMRX 
IFCVXMIN.NE.O. .OR. VXMRX. NE.O.) GO TO 30 
VXMIN-VPNOMC 1 , KNOM) 

VXMRX- VPNOM C 2 , KNOM ) 

VYMIN-VPN0MC3, KNOM) 

VYMRX-VPN0MC4, KNOM') 

30 CONTINUE 

XMIN-O. 

XMAX-O. 

WRI TEC 1 ,401 ) 

READC 1 , # ) XMIN, XMRX, YMIN, YMAX 

IF (XMRX.NE .0. .OR. XMIN. NE.O.) GO TO 31 

XMIN-O, 

XMAX-500. 

YMIN--1 00. 

YMAX-900. 

31 CONTINUE 

XLOOK-O. 

WRI TEC 1,404) 

READ Cl,#) XLOOK 
IF (XLOOK. EQ.O. ) XLOOK- 1 . 


NCRV-4 

WRI TEC 1,402) NPLT 
RERD Cl,#) NCRV 
IFCNCRV.EQ.O) GO TO 32 
C 

C NCRV - 0 IS AN ABORT 
C 
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0131 

0132 

0133 

0134 
0135- 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0148 
0147 
0140 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 
0150 

0159 

0160 
0161 
01 62 

0163 

0164 

0165 

0166 
0167 
0160 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 
0100 
0101 
0182 
0103 
010 4 
0165 
01 86 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 


CRLIB-0. 

WRITE (1,403) 

RERD( 1 , #) CRLIB 
IFCCRLIB.EQ.O. ) CRLIB". 2 


32 


C 


CONTINUE 
MCRV "NCR V 


IFRRf!-1 
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DO S3 ICRV»1,MCRV 
IF(NCRV.EO.O) GO TO 44 
NRVG-NRVGDF 
WRITE (1,501) NRVGDF 
READ(1,#) NRVG 

I F ( NRVG . E 0 . 0 . AND . I CRV . E Q . MCRV . RND . NPLT . E Q . NPLO T ) NCR V -0 
IF (NRVG. EQ. 0. AND* ICRV. NE . MCRV ) GO TO S3 

IF (NRVG. EQ. Ot AND. I CRV. EG. MCRV, RND* NPLT. NE . NPLOT) GO TO 93 
IF(NCRV.EQ.O) GO TO 44 


C 

C NRVG-0 IS R 4# CURVE ## ABORT 
C 

IF (LU. EQ. 32 ) WRI TE( 1 , 502 ) 

MTQPT-0 

IRNS-1HN 

WRITE (1,504) 

READ (1,505) I RNS 
IF (IANS. NE . I YES) GO TO 33 


C 


MTOPT-1 
MTCNT "MTCNT+1 
WRI TE( 1,510) 

CRLL LXEC(0,PGSTP) 

105 CONTINUE 
IRNS-1HN 
I0PEN-2H0P 
WRITE (1,506) 

READ (1,505) I RNS 
IF (IANS.EQ.IYES) IQPEN-2HBL 
33 CONTINUE 


IFdFRRM. NE. 1 ) GO TO 40 
CALL EXEC ( 1 1 , ITIME, IYER) 

IRM-2HRM 

I F ( I TI ME ( 4 ) . GE . 1 2 ) IAM-2HPM 
IFCITIMEC4). GT .12) ITIME (4 ) «I TIME (4 ) - 1 2 
I HOUR "ITIME (4 ) 

MI NUT" I TI ME ( 3 ) 


IF(MQD(IYER,4) .EQ.O) IDAYS(2)-29 
I DUMB -0 
DO 35 1-1,12 
M7-I 

I DUMB ■ I DAYS (I ) +1 DUMB 
I F ( I TI ME ( 5 ) . LE . I DUMB ) GO TO 36 

35 CONTINUE 

36 MON»M? 

I DAY" I DAYS (M7 ) - (I DUMB -I TIME ( 5 ) ) 

I YEAR"I YER-1 900 

IFdFRRM. NE. 1 ) GO TO 40 
CALL EXEC(8, PTLID) 

103 CONTINUE 
IFRRM-0 
WRITE (1,507) 

READ( 1 , #) PAUSE 

40 CONTINUE 13 


c 


0197 

0190 

0199 

0200 
0201 
0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 
021 1 
0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 
0221 
0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 


C 


C 


C 

C 


c 

c 


CALL EXECC8, BISEG) 

101 CONTINUE 

IF (MTOPT.EQ. 1 ) CALL EXECC8, TRASH) 
104 CONTINUE 


44 CONTINUE 

CALL EXECC0,PTLID) 

102 CONTINUE 

IFCNCRV.EQ.O) GO TO 41 
93 CONTINUE 

IFCNCRV.EQ.C) GO TO 41 
92 CONTINUE 
GO TO 91 
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90 CONTINUE 
WRI TE ( 1 993 ) 

990 FORMAT (1 X, 14, " FILES WRITTEN TO TAPE") 

999 FORMRTC IX, "IMPROPER LU OR ID... USE PARAMETERS IN RU COMMAND") 

91 CONTINUE 


41 CONTINUE 

WRI TEC 1,990) MTCNT 
END 

BLOCK DATA XYZ 

COMMON/XNY/X C 5 1 2 ) , Y C 5 1 2 ) , I TOUT, NRVG, IFRRM, I OPEN, KFILE 
COMMON/GRRPIVI GCB < 1 92 ) , I PRAM ( 5 ) , I OBUF C 20 ) , I DCB ( 1 44 ) , I FI LE C 3 ) 
XMIN, XMAX, YMIN, YMAX, NCRV, ICRV, VXMIN, VXMRX, VYMIN, VYMflX, NPLT 
- , NFLCT, CALIB, XLOOK, I ENTRY, MTOPT, MON, JURY, I YEAR, I HOUR, MI NUT, I AM 
COMMON I NSTN1 , I NSTN2 , I NSTN3 , I NSTN4 , I NSTN5 


END 

END# 
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0021 
0022 
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0024 

0025 

0026 
0027- 
0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 
003? 

0038 

0039 

0040 

004 ; 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 


r TN4,L 

PROGRAM BISEG(5 > 

CQMMON/'XNY/X ( 5 12), DRTR( 5 1 2 ) , I TOUT, NflVG, IFRAM, I OPEN, KFILE 
CQMMON/GRRPH/1 GCB C 1 92 ) , I PRRM ( 5 > , 1 0 BUF < 2 0 ) , X DCB ( 1 4 O , I F I LE ( 3 > 
-,XMIN„ XMRX , YMIN, YMRX , NCRV , ICRV, VXMI N , VXMRX , VYMI N , VYMRX , NPLT 
-, NPL3T, CRLIB, XLOOK, IENTRY, MTOPT, MON, I DRY, I YERR, I HOUR, MI NUT, IRM 
COMMON I NSTN1 , I NSTN2 , I NSTN3 , I NSTN4 , I NSTN5 
DIMENSION I DflTR < 1 02 6 > , I RE G ( 2 ) , I CBUF ( 5 ) , I BUFF ( 1 5 ) 

DIMENSION I WRTB ( 2 ) , ZERO (512), STORE (512), XDRTR (512) 

DIMENSION IDRT2C514) 

EQUIVALENCE (REG, IREG( 1 ) ) 

C 


I BUFF < 1 ) “MON 
I BUFF (2) -I DRY 
I BUFF (3 )*I YEAR 
IBUFF (4 ) "IHOUR 
I BUFF (5 ) "MI NUT 
IBUFF (6) -IRM 
IBUFF <7 ) -NflVG 
IBUFF (8 ) -XLOOK 
I BUFF ( 9 ) ■ 1 0#CALI B 
IBUFF ( 1 0) -I OPEN 
IBUFF (1 1 ) -KFILE 
IBUFF < 1 2 ) -0 
IBUFF < 1 3 ) -0 
IBUFF < 1 4 ) -0 
IBUFF ( 1 5 ) -0 
C 
C 
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THIS SEGMENT READS DRTR FROM THE BIQMRTION. IT AVERAGES 
NflVG SHOTS, ERCH CONSISTING OF THE FIRST 512 POINTS OF 
BIOMRTION MEMORY. THE SECOND 512 POINTS OF BIOMRTION 
MEMORY ARE USED TO DETERMINE R BASELINE (I.E. ZERO) 

LEVEL. THE LAST 1024 POINTS RRE DISGRRDED TO ENRBLE 
RERDING DflTR RT R 2 HERTZ REP RATE. TRYING TO READ 
RLL 2048 WORDS SLOWS THE TRANSFER TIME TOO MUCH RND 
R 1 HERTZ MAXIMUM IS ABOUT RLL YOU. CRN GET FROM THIS 
SEGMENT. THIS SEGMENT USES THE DVM72 DRIVER SUPPLIED 
BY HP FOR THE 12566B 1SF . REFER TO THE DRIVER MANUAL 
FOR DETRILS CONCERNING THE VARIOUS EXEC CALLS. 

THIS SEGMENT IS BASED ON THE STAND-ALONE PROGRAM " YBIOM" 
WHICH IS USED SIMPLY TO READ DATA AND PRINT OUT THE AVERAGE. 


NCRV-0 IS AN ABORT 


IF (NCRV. EQ. 0) GO TO 41 
C 


N1 026 -1026 


C 

C N IS THE NO. OF DATA WORDS + NO . OF WORDS IN COMMAND STRING 
C 

ICWTO-1 33B 
C 

C ICWTO IS THE CONTROL WORD FOR SETTING TIME OUT 
C BIT 6 SAYS RESETTING TIME OUT VALUE 
C 

IF (I TOUT . L£ . 20) ITOUT-50 
C 

C I TOUT -50 SAYS T.O. VALUE IS TO BE 500 MILLISECONDS 
C 

CALL EXECC3, ICWTO, ITOUT) 

C 15 
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0071 
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0087 
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0096 
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0098 
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0100 
0101 
0102 
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0108 

0109 
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0111 
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0115 

0116 
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0123 

0124 

0125 

0126 
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0130 


C 

C 

C 

C 

C 

C 

C 

C 

c 


IDATA<1 )-14 
ICBL-2 
I CBUF < 1 ) -5 
I CBUF (2 ) "6 
ICODE-1 
ICW2-1 1533B 

15 (14 
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QUALITY 


13 


1 


12 ) 

1 


<11 10 


1 


S) 

1 


(8 

1 


6 ) 

1 


(5 


3) 

1 


(2 


1 
1 
3 

WITH 


0) 

1 


ICW2 IS THE CONTROL WORD FOR °ERFQRMING COMMAND SEQ. 

BIT 12 SAYS USING I CBUF 

BIT 10 SAYS T. 0. IS LEGAL. MUST BE SET FOR IDATA<1)»12 
BIT 9 SAYS SUBFUNCTION IS IN I DATA ( 1 ) 

BITS 8 AND 6 SHY STANDARD READ. 


INTERUPT 


IWRTL-2 
I WRTB ( 1 ) » 1 
IWRTB<2)-0 
ICWRT-2333B 

C 15 <14 13 12) <11 10 9) <6 76) <54 3) <2 1 0) 

C 1 111111 

C 2 3 3 3 

C I CWRT IS THE CONTROL WORD FOR UNINTERUPTED WRITE 
C BIT 10 SAYS T. 0. IS LEGAL 

C BITS 7 RND 6 SRY WRITE WITHOUT INTERUPT UNTIL COMPLETE 
C THE TWO DATA WORDS IN I WRTB EFFECTIVELY ISSUE AN OSR 
C 

C ■ *************************** 

C 

C RESET AFTER EACH DATA TRANSFER USING I CBUF IN READ EXEC CALL 


C COMMAND 5 IS A CLF 
C COMMAND 6 IS AN STF 

C SUBFUNCTION 14 RLLOWS CLC, STF IN FINI 
C 


101 F0RMAT<1X,I4,3X,8F6.1,4X,06,F7.1 ) 
DO 30 1-1,512 
STORE <I) -0.0 
30 CONTINUE 

C ******** , 

L-1 


C ******** 

ICTRP-10B 

IF (MTQPT.EQ. 1 ) CALL EXEC(2, ICTAP, IBUFF, 1 5 ) 

C 

20 CONTINUE 

REG-EXEC < I CODE , I CW2 , 1 DATA, N1 026 , 1 CBUF , I CBL ) 
I AREG-IREGC 1 ) 

REG- EXECC2, ICWRT, IWRTB, I WRTL) 

C ******* 

DO 22 JHU-1,512 
XDATA< JHU)-IDATA( JHU+2 ) 

22 CONTINUE 
C ******* 


IDRT2 < 1 ) -L 
IDRT2 <2 ) -1000 
DO 10 1-1,512 
DATA < I ) - 1 DATA ( I “i*2 ) 

ZERO < I ) -I DATA <1+51 4 ) 

10 CONTINUE 

CALL VSUM(ZBRR, ZERO (1), 1,51 2) 

ZBAR--1 .0*ZBRR/512. 

CALL VS AD < ZBAR , DATA , 1 , DATA ,1,512) 

CALL VADD(DATA, 1 , STORE, 1 , STORE, 1,512) 

WRITEU, 101 ) L, (DATH<I ) , I -100, 450,50) , IRREG, ZBAR 
C ******** 

IF (IAREG. EQ, 35004B) GO TO 20 
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0131 

0132 

0133 

01 34 

0135 
01.36 
0137 
0130 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 
01 40 
w’l 49 

0150 

0151 

0152 

0153 

0154 


IF (MTQPT.NE. 1 ) GO TO 32 
DO 33 1-3,514 
I DRT2 < I ) - 1 DRTR <1 ) +ZBFIR 
33 CONTINUE 

CRL.L EXECC2, ICTRP, IDRT2 ,514) 

32 CONTINUE 
L-L+1 

IFCL.LE.NRVG) GO TO 20 
21 CONTINUE 
ICEOF-1 10B 

IF CMTQPT.EQ. 1 ) CALL EXEC(3, ICEOF ) 
C ###*#### 


C 

C 


40 

41 
100 


XRVG-1 .0/FL0RT<NRVG> 

CRLL VSMYCXAVG, STORE, 1 , DRTR, 1,512) 

DO 40 IJK-1,490, 10 

IMIN-IJK 

IMRX-IMIN+9 

WRITEC1, 100) (DRTR (I ),I-IMIN, IMRX) 

CONTINUE 
CONTINUE 
FORMAT (10F7.0) 

GO TO INSTN1 
END 

END * ORIGINAL PACE'S 
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0026 

0027 

0028 
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FTN4 # L 

PROGRAM TRASH (5 ) 

COMMON /XNY/X C512),Y<512),I TOUT, NAVG, I FRAM, I OPEN, KFILE 
COMMON I NSTN1 , 1 NSTN2 , 1 NSTN3 , 1 NSTN4 , 1 NSTN5 
DIMENSION I BUFFI < 15),IBUFF2<514) 

C 

C DATE 

I BUFFI ( 1 ) "13 
I BUFFI C2J-39 
I BUFFI (3 ) *83 


C 

C TIME 

I BUFFI (4)-13 
I BUFFI <5)-71 
I BUFFI (6 ) *2HXX 
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C 

C NUMBER OF SHOTS 
I BUFFI (7)-1 


C 

C HOROZONTAL/VERTICAL 
I BUFFI (8 )“Q 


C 

C TIME /POINT IN TENTHS OF MICROSECONDS 
I BUFFI <9 ) "0 
C 

C OPEN/BLOCKED 

IBUFF1 < 10)»2HXX 

IBUFF1 (11) "KFILE+1 
t oi ice i 1 1 o \ s n 

o> t/MM I I % • &» * ”V 

I BUFFI < 1 3) "0 
IBUFF1 <14)-0 
I BUFFI ( 1 5) “0 
IBUFF2C1 )-1 
IBUFF2 (2) *999 
DO 10 I “3, 514 
IBUFF2CI ) *0 
10 CONTINUE 
ICNWD-10B 

CALL EXEC(2,ICNWD,IBUFF1 ,15) 

CALL EXECC2, ICNWD, IBUFF2, 514) 

ICNWD-110B 

CALL EXECO, ICNWD? 

ICBK-1410B 

ICFWD-1310B 

CALL EXECC3, ICBK) 

CALL EXECC3, ICBK) 

CALL EXECC3,ICFWD) 

GO TO I NSTN4 ‘ 

END 

END! 
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0001 

FTN4, 

L 


0002 



PROGRRM PQ5TPC5) 


0003 



COMMON /XNY/X (512),Y(312),I TOUT, NRVG, I FRRM, I OPEN, KFILE 

0004 



COMMON I NSTN1 , I NSTN2 , I NSTN3 , I NSTN4 , I NSTN5 


0005 



DIMENSION I BUFFI ( 1 5 ) , I BUFF2 (51 4 ) 


ooos 

C 




000? 



ICNWD-10B 


0008 



ICFWD*1 31 OB 

ORfGffWHL PAQg Jg 

0009 



ICBK-1410B 

0010 



IREWND-410B 

0F POOR QUALITY 

0C1 1 

C 




0012 



WRITE ( 1,200) 


0013 


200 

F0RMAT(26X," SHOTS H/V T/P OP/BL” ) 


0014 

C 


CALL £XEC(3, IREWND) 


0015 



KOUNT-O 


0016 


95 

CALL EXECC1 , I CNWD, IBUFF 1 ,15) 


001? 



KFILE -I BUFFI (11 ) 


0010 



KOUNT-KOUNT-M 


0019 



WRITE (1,1 01 ) KOUNT, (IBUFF 1 (I), I *1,15) 


0020 


101 

FORMAT ( 2X, 21 4 , " /“ , I 2 , " ' ✓« , 1 2 , 2X, 1 2 , « S " , 1 2 , 

1 X, A2, 315, 3X, A2, 515 

0021 

C 




0022 

C 

CHECK TO SEE IF TRRSH FILE. I.E., MONTH RND 

HOUR ARE 13. 

0023 

C 




0024 



IFdBUFFI ( 1 ) ,EQ. 1 3 .AND. IBUFF1 (4) ,EQ. 13) 

GO TO 104 

0025 



C. .1 EXEC(3, ICFWD) 


0026 



GO TO 95 


0027 


104 

CONTINUE 


0028 



CALL EXECO, ICBK) 


0029 



CALL EXEC (3, ICFWD) 


0030 



GO TO INSTN5 


0031 



END 


0032 



ENDt 



j 
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C 

C 

C 


PROGRAM PTLID(5> 

COMMQN/XNY/'X ( 5 1 2 ) , Y < 5 1 2 ) , I TOUT , NflVG , I F RAM , X OPEN , KFI LE 
COMMQN/GRRPH/I GCB ( 1 9 2 ) , I PRAM (5 ) , IQBUF (20) , IDCB( 1 44 ) , IFILE (3 ) 

- , XMX N , XMRX , YMIN, YMRX , NCR V , I CRV , VXMI N , VXMRX , VYMI N , VYMAX , NPLT 
-,NPLOT, CALIB, XLQQK, I ENTRY, MTOPT , MON, I DAY, I YEAR, I HOUR, MI NUT, I AM 
COMMON INSTN1 , INSTN2, I NSTN3 , I NSTN4 , I NSTN5 
DIMENSION DISK 6) 

100 FORMAT C SPO" ) 

101 FORMATS SP1 " ) 

102 FORMAT <" SP2" ) 

103 FORMAT (" SP3" ) 

104 FORMAT (" SP4” ) 

105 FORMAT C " SP5" ) 

106 FORMAT <" SP6" ) 

107 FORMAT ("SP?" ) 

108 FORMAT ("SP8" ) 

IENTRY-IENTRY+1 
NP-512 

DO 10 1-1, NP 
XCD-I 
10 CONTINUE 
LU-IPRAM( 1 ) 

ID-IPRAM(2 ) 

IF (ID. EQ, 1 .AND. I ENTRY. EQ. 1 ) 

-CALL PLQTRdGCB, ID, 4, LU, IOBUF, 20) 

IF (ID.EQ.2. AND.IENTRY.EQ. 1 ) 

-CALL PLO TR ( I GCB , I D, 4 , LU, 1 0 BUF , 20 ) 

NCRV.FQ ; Q IS AN ABORT 

IF(NCRV.EQ.O) GO TO 41 


C 

C 

c 

c 

c 

c 

c 


IP-MOD<ICRV, 8 ) 

IFCIP.LE.O) IP-0 

IF(IP .EQ. 1 . AND. LU. NE . 32 ) WRITE (LU, 1 01 ) 
IF(IP .EQ. 1) WRITE (LU, 101) 

IF(IP .EQ. 2) WRITE(LU, 102) 

IF (IP .EQ. 3) WRITE(LU, 103) 

IFdP .EQ. 4) WRITE(LU, 104) 

IFdP .EQ. 5) WRITE(LU, 105) 

IFdP .EQ. 6) WRITE(LU, 106) 

IF(IP .EQ. 7) WRITE(LU, 107) 

IFdP .EQ. 8) WRITECLU, 108) . 

IFdP .EQ. 2 . AND. LU. NE . 32 ) WRITE (LU, 1 02 ) 
IF(IP .EQ. 3 . AND . LU . NE . 32 ) WRITE (LU, 1 03 ) 
IFdP .EQ. 4 . AND. LU. NE . 32 ) 

IF(IP .EQ. 5 . AND. LU; NE . 32 ) 

IFdP .EQ. 6 . AND. LU. NE . 32 ) 

IFdP .EQ. 7 . AND. LU. NE . 32 ) 

IFdP .EQ. 9 . AND. LU. NE . 32 ) 

CALL SETAR ( I GCB ,1.5) 

CALL VI E WP d GCB, VXMI N , VXMHX , VYMI N, VYMAX ) 

CALL WI NDW (I GCB , XMI N , XMRX , YMI N , YMRX ) 

HGT-3.0 

AR-O. 4 

SANG -0.0 

ICOOR-O 

CALL CSI 2E ( I GCB, HGT, AR, SANG, I CO OR 1 ) 

CALL FXDdGCB, 1 ) 

XTIC- (XMIN-XMRX)/'S . 

YTIC-(YMAX-YMIN)/5. 

IF (IFRRM.NE. 1 ) GO TO 42 
CALL FRAME (I GCB) 


WRITE (LU, 104 ) 
WRITE (LU, 105) 
WRITE (LU, 106) 
WRITE (LU, 10?) 
WRITE(LU, 108) 
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0065 

0066 
0067 
0066 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 
0076 
0079 
0060 
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0062 

0083 

0084 

0085 

0086 
0067 
0088 

0089 

0090 

0091 

0092 

0093 

0094 
0035 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
01 1 1 
0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 
0127 


C 


C 

C 


CALL FRAME Cl GCB) 

CALL LAXESCIGCB,XTIC, YTIC, XMIN, YMIN, 1 . , 1 . ,2. ) 
VXMI3-VXMIN-3. 

VXMA3 • VXMRX+6 . 

VYMI3-VYMIN-3. 

VYMA3 -VYMRX+3 . 

CALL VI E WP (I GCB , VXMI 3 # VXMA3 , VYMI 3 , VYMA3 ) 

CALL WI NDW C I GCB , VXMI 3 , VXMR3 , VYMI 3 , V YM A3 ) 

DO 43 IDIST-1,6 

DIST(I DIST) ■ CXMIN- < I DI ST-1 )#XTIC)#CALIB*. 3/2 
43 CONTINUE 

CALL LORGCIGCB, 1 ) 

CALL MO VE (I GCB , VXMI 3 , VYMAX ) 

CRLL LABEL (I GCB) 

WRITE CLU, 431 ) CDISTCIDIST) ,IDIST»1 ,6 ) 

431 FORMAT C 1 X, 5 CF4 . 1 , 6X) ,F4 . 1 ) 
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VXMA2-VXMR3-3. 

CALL MOVE CIGCB,VXMR2, VYMI N) 

CRLL LDIRCIGCB, 1.57) 

IF (XLOOK. EQ. 2 ) GO TO 30 
CRLL LABEL (I GCB) 

WRITE CLU, 432) NAVG, MON, I DAY , I YEAR, I HOUR , MI NUT, I AM 
GO TO 31 

30 CONTINUE 

CALL LABEL (I GCB) 

WRITE CLU, 433 ) NAVG, MON, I DRY, I YEAR, I HOUR, MINUT, I AM 

31 CONTINUE 

CRLL LDIRdGCB, 0. ) 

432 FORMAT d 4," SHOT AVERAGE - HORIZONTAL" ,x, 

-6X,I2," /" »r2."/" ,I2,4X,I2," *" ,12, 1X,R2> 

433 FORMAT (I 4, " SHOT AVERAGE - VERTICRL 

-6X,I2," s%12," , 12, 4X, 12, " : " , 12, 1 X, A2 ) 

CRLL MO VE d GCB , VXMI N , VYMI N ) 

CALL VIEWPdGCB, VXMIN, VXMAX, VYMIN, VYMRX) 

CRLL WINDWdGCB, XMIN, XMAX, YMIN, YMAX) 


IFCIFRAM.EQ.1 ) GO TO INSTN3 , 

42 CONTINUE 

DO 20 J " 1 , NP 

IFCJ.GT.1 ) GO TO 21 

CALL MOVE (IGCB, X( J ) , YC J ) ) 

GO TO 20 

21 CALL DRAW C I GCB , X C J ) , Y C J ) )' 

20 CONTINUE . 

XTI ME - XMI N+3 . 2 * ( XMAX* xfll N)/5 . 
YTIME-YMAX-ICRV*CYMRX-YMIN)/10. 

CALL MOVE Cl GCB, XTI ME, YTIME ) 

rai i i mnri ftcpni 

IFCMTOPT.NE. 1 ) WRITE CLU, 201 ) IHOUR, MINUT, I AM 
IF CMTOPT.EQ. 1 ) WRITE CLU, 202 ) IHOUR, MINUT, I AM, KFILE 

201 FORMATC IX, 12, " S" ,12, 1X,A2) 

202 FORMATC IX, 12," ,12, 1X,A2, IX, 14) 

CRLL XMITCIGCB) 

41 CONTINUE 

IF CNPLT. EQ. NPLOT. AND. ICRV. EQ.NCRV) WRITE CLU, 1 00) 
IFCNCRV.EQ. 0) WRITE CLU, 1 00) 

IFCNPLT.EQ. NPLOT. AND. ICRV. EQ.NCRV) CALL PLOTRCIGCB, ID, 0) 
IFCNCRV.EQ. 0) CALL PLOTRCIGCB, ID, 0) 

GO TO INSTN2 

END 

END? 
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0022 
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0024 

0025 

0026 

0027 

0028 
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PROGRAM TEMPE 

COMMON/XNY ■'NR VG , I F RAM 1 1 BUF 1 ( 1 5 > , I BUF 2 ( 5 1 4 ) 

COMMON I YES, NO 

COMMON/I TX C5/I TIC1 , ITIC2, 1 TICS, ITIC4, IT2C5, ITI C6 
COMMON IREWND,ICNWD,ICNWD1 , ICNWD2, ICNWD3, NPTS 
COMMON/GRAPH/I GCB ( 1 92 ) , IPRRM(5 ) , IQBUF (20) , I DCB( 1 44 ) , IFILE (3 ) , 
-XMIN, XMRX, YMIN, YMRX, NCRV, ICRV, VXMIN, VXMAX, VYMIN, VYMAX, NPLT, 

- NPLO T, CALI B, I ENTRY, I RM 
COMMON I NSTN2 , I N9TN3 , I NSTN4 , I NSTN5 
DIMENSION VPN0M(4, 4 ) , XDRYSC 1 2 ) , ITIME (5 ) 

INTEGER PTTEM(3 ) , FNTEMC3 ) 

EQUI VRLENCE (IPRRM( 1 ),LU), (IPRRM(2 ), ID) 

DRTR VPNOM/1 0. , 70. , 57 . , 97. , 1 0. , 70. , 7 . , 47 . , 85 . , 1 45 . , 57 . , 97 . , 
-85. ,145. ,7. ,47./ 

DRTR IDAYS/31 ,28,31 ,30,31,30,31,31,30,31 ,30,31/ 

IREWND-410B 

ICNWD-10B 

ICNWD1-1310B 

ICNWD2-1410B 

I CNWD3 * 1 1 09 

NPTS-514 

ITIC1-1H: 

ITIC2-1H- 

ITIC3-1HH 

ITIC4-1HV 

ITICS-2H0P 

ITIC6-2H&L 

IYES-1HY 

N0-1HN 


RSSIGN 

ASSIGN 

ASSIGN 

RSSIGN 


TO 

TO 

TO 

TO 


INSTN2 
INSTN3 
INSTN4 
INSTN5 
THE 


LU AND ID OF THE PLOTTING DEVICE 


COMBINATION IS USED 


68 

42 

600 


601 


602 


603 


351 


102 

103 

104 

105 

THIS CALL PICKS UP 
CALL. RMPARCIPRRM) 

THIS ASSURES THAT A PROPER LU/XD 
I F C LU . LE . 0 . OR - 1 D . LE . 0 ) GO TO 90 
IF'ID.EQ.I )LU-1 

IF(ID.EQ.2. AND. LU . NE . 32 )GQ TO 90 
IFCID. GT. 2 )G0 TO 90 
WRITEd ,600) 

FORMAT (R1 ) 

FORMAT (" YOU MUST ANSWER Y(YES) OR N(NO) ONLY II, TRY AGAIN" ) 
FORMATC" THIS PROG. IS THE MAIN PROG. FOR SETTING UP PLOTS",/, 

" IT READS IN THE PLOT DEVICE LU AND ID THRU. RMPAR. . . 

" I. E. RU, STATS, 1,1 PLOTS ON THE CRT OR",/, 

" RU,STRTS,32,2 PLOTS ON THE PLOTTER",/) 

WRITE (1,601) . 

FORMATC QUESTIONS WILL BE ASKED CONCERNING ' HOW MANY PLOTS, 

" SCALING OF EACH INDIVIDUAL PLOT AND PEN NUMBER" ) 

WRITE (1,602) 

FORMATC THIS MAIN PROG. ASSUMES THAT THE ACTUAL PLOTTING",/, 

" ROUTINE AND TAPE POSITIONING ARE PROGRAM SEGMENTS",/) 
WRITE ('! ,603) 

FORMATC THE MAIN 
" AND THEN 
" ON TAPE, 

" YOU MUST 
" ON, THRU. 

CONTINUE 
FNTEMd )-2HFN 
FNTEM(2)-2HTE 

FNTEM(3)-2HM 22 




REWIND THE MAG. TAPE" , /, 
ALL THE DATA FILES " 


PR0<3. STATS WILL FIRST 
PRINT OUT HEADINGS OF 
UP TO A TRASH FILE 
INPUT THE FILES OF THE FILE YOU WISH 
THE KEYBOARD TO ACCESS THAT FILE.") 






TO WORK" 






? 


0065 

0066 

0067 

0068 
0063 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0081 
0082 

0083 

0084 

0085 

0086 
0087 
0080 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
01 1 1 
0112 

0113 

0114 

0115 

0116 

0117 

0118 
0113 
0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 


C 

C 


C 

C 


QRfQifWU, P&Qg fsr 
p 00R QUALITY 


CALL EXEC(8,FNTEM) 

104 CONTINUE 

CALL EXECd,ICNWD,IBUF1 , 15) 

NRVG-IBUF1 (7) 

CRLI B •FLOAT ( 1 BUF 1 (9))/10, 

WRITE (1,273) 

273 FORMRT(/, " THE FOLLOWING IS HERDING OF YOUR DRTR FILE",/) 
WRITE (1,274) 

274 FQRMRK0X, " 4 SHOT" , 4X, " DRTE" , 6X, " TIME" , 6X, " H/V" , 7X, 

T/P" , 7X,"QP/BL" ,/,66(" •" ),/) 


~l* 


IF (I BUF 1 (0).EQ.1 )IBUF1 (8) -I TICS 
IF (I BUF 1 (0 ) . EQ . 2 )I BUF 1 (8)-ITIC4 
IF (I BUF 1 (10) .EQ. 1 )IBUF1 (10) -I TICS 
IF(IBUF1 (10).EQ.2)IBUF1 (10)-ITIC6 

WRI TE ( 1 , 275 ) I BUF 1 (7),IBUF1 (1 ),ITIC2,IBUF1 C2),ITIC2,IBUF1 (3), 
-IBUF1 (4),ITIC1 , IBUF1 (5),IBUF1 (6),IBUF1 (0),IBUF1 (9), IBUF1 ( 10) 
275 F0RMRT(1OX,I4,3X,I2,R1,I2,A1,I2,3X,I2,R1 ,I2,1X,R2, 
-4X,R1,0X,I3,9X,R2) 


WRITE (1,507) 
READ (1,4) PAUSE 


DO 20 J-1,NRVG 

CRLL EXECd ,ICNWD,IBUF2,514) 

DO 60 1-1,514 

BUF2 ( I ) -FLOAT < I BUF 2 (I ) ) 

60 CONTINUE 

CRLL VRDDCBUFSU, 1 , BUF2, 1 , BUFSU, 1,514) 

« A Aftklf Y Lll If* 

/ Vi LLini ■ t mLIP 


CRLL £XEC(3, I CNWD1 ) 

CRLL VSMY(1 ./NRVG, BUFSU, 1, BUFSU, 1,51 4) 

I ENTRY -0 
NPLOT-O 
WRITE (1,300) 

300 FORMRTC HOW MRNY FRAMES/SCREEN DO YOU WISH TO HAVE? <DEFAULT *4>" ) 
RERD( 1 , #)NPLOT 
IF(NPLOT.£Q.O)NPLOT-4 
PTTEM( 1 )-2HPT 
PTTEM(2 ) -2HTE 
PTTEM(3 ) -2HM 
DO 92 NPLT-1 , NPLOT 
IF (NPLT . EO. 1 )G0 TO 401 
020 WRI TE (1 , 1 06 ) 

106 FORMRTC WISH TO RNRLY2E RNOTHER PORTION OF SAME FILE (Y/N)?" ) 
READC1 ,60)IRNS 
IT <IANS.EQ.IY£S)GQ TO 401 
IF (IANS. SQ. NO )G0 TO 310 

IF( (IRNS.NE.IYES) . OR. (IANS. NE .NO ) ) WRITE d , 42 ) 

GO TO 020 
401 VXMIN-O. 

VXMRX-O. 

KNQM-MQDCNPLT, 4 ) 

IF (KNQM.EQ.O)KNOM-4 

WRI TEd ,400) NPLT, ( VPNOMC JNOM, KNOM) , JNOM-1 , 4 ) 

400 FORMRTC INPUT VIEWPORT COORD. : VXMIN, VXH AX, VYMIN, VYMAX" , /, 

-"•FOR FRAME ",I3," <DEFAULT VRLUES-" , 3 <F3 . 0, " , " ) , ^3 . 0, " >" ) 

RERDd , *) VXMIN, VXMRX, VYMIN, VYMAX 
r c ( VXMI N . NE . 0 . . OR . VXMRX . NE . 0 . ) GO TO 30 
vWlIN-VPNQMd ,KNQM) 

VXMRX- VPNOM( 2; KNOM) 

VYMIN-VPNQM(3 , KNOM) 

VYMAX- VPNOMC 4, KNOM) 

30 CONTINUE 
XMIN-O. 


23 


0131 

0132 

0133 

0134 

0135 

0136 
013? 

0138 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 
014? 
0148 
0145 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 
0150 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 
0188 

0187 

0188 

0189 

0190 
0131 

0192 

0193 

0194 

0195 

0196 


XMRX-500. 

YMIN-- 1 00, 

YMRX-900, ORIGINAL PAGE IS 

31 CONTINUE OF POOR QUALITV 

RLTMX- ( XMI N-5# C ( XMI N-XMRX > /5 . ) ) #CALI B# , 3/2 y MU 1 Y 

38 WRI TE < 1 , 4 1 ) RLTMX 

41 FORMRTC/," THE MAX. ALT . FOR THIS FILE IS-" ,F5 .2," KM.",/) 

WRITEC1 ,22) 

22 FORMRTC" WHPT RRNGE OF RLT. YOU WISH TO LOOK RT <IN KM.> I.E", 

< 1 .5,3.5)" ) 

RERDC 1 , #)RLTM1 , RLTM2 

I F < < RLTM1 , GE . 0 . ) . OR . < RLTM2 . LE . RLTMX ) ) GO TO 45 
WRITE <1,44) 

44 FORMRTC/," RLT. RRNGE ERROR il, TRY RGRIN" ) 

GO TO 38 

45 NRLTM1 -IFIX<FLQAT<NPTS-2)#<RLTM1 /RLTMX) ) 

NI ILTM2 -I FI X < FLORT < NPTS-2 ) # < RLTM2/RLTMX ) > 

XMIN- CFL0RTCNRLTM1 )/FLQRTCNPTS-2 ) )#S00. 

XMRX- CFL0RTCNALTM2) /FLOAT CNPTS-2)) #500. 

WRI TE < 1 , 1 0? ) NRLTM1 , NRLTM2 , XMI N, XMRX, CRLI B 

10? F0RMRT(2X, 15, 5X,I5,5X, F6.0, 5X, F6 . 0, 5X, FB . 3 ) 

WRITE (1,50?) 

RERD Cl-,#) PAUSE 
NCRV-2 
MCRV-NCRV 
DO 93 ICRV-1 , MCRV 
WRI TEC 1,776) ICRV,MCRV 
776 FORMRTC" ICRV-" ,15," MCRV-" ,15) 

I F Cl CR V . EQ . 2 ) CALL SMO 0 T ( BUF SU , NPTS , I CR V ) 

IF CLU. EQ. 32 )WRITE C 1 , 507 ) 

50? FORMRTC" ## PRESS RETURN TO CONTINUE ##" ) 

RERD Cl,#) PAUSE 
IFCICRV.NE. 1 )G0 TO 77 
CALL EXECC 1 1 , I 71 ME , I YER) 

IRM-2HRM 

IF CITIME ( 4 ) . GE c 1 2 )IRM-2HPM 
IFCITIMEC4). GT. 1 2 )ITIME <4 ) “ITIMEC4 ) -1 2 
IH0UR-ITIMEC4) 

MI NUT -ITIME C3) 

IF CMODCIYER, 4 ) .EQ.0)IDRYS(2)-29 
IDUMS-0 
DO 397 1*1,12 
M7-I 

I DUMB -I DAYS (I )+IDUMB 
I F C I TI ME ( 5 ) . LE . I DUMB ) GO TO 378 
397 CONTINUE 
378 M0N*»M7 

IDRY-IDRYSCM7 ) - CIDUMB-ITIME C5 ) ) 

I YERR-I YER-1 900 
IFCICRV.NE. 1 )G0 TO 77 
IFRRM-1 

CALL EXECCB,PTTEM) 

103 CONTINUE 
IFRRM-0 

77 IFCICRV. EQ. 1 )GQ TO 93 
CALL EXEC(8,PTTEM) 

102 CONTINUE 
93 CONTINUE 
92 CONTINUE 

318 WRI TE C 1 ,319) 

319 FORMRTC/," DO YOU WISH TO ANALYZE ANOTHER SET CY'h)?" ) 

RERDC1 ,S8)IRNS 
IFCIRNS. EQ. I YES)GO TO 351 
IF (IANS. EQ. NO )G0 TO 91 

IFC ClflNS.NE.IYESJ.flND. CIRNS. NE . NO ) ) WRITE C 1 , 42 ) 

GO TO 318 


24 


019? 

0198 

0199 

0200 
0201 
0202 

0203 

0204 

0205 
0208 

0207 

0208 

0209 

0210 
0211 
0212 

0213 

0214 

0215 

0216 
021 ? 
0218 

0219 

0220 
0221 
0222 

0223 

0224 

0225 

0226 


90 CONTINUE 
WRITE (1,999) 

999 FQRMRTC" IMPROPER LU OR ID.... USE PRRRMETERS IN RU CQMMRND" ) 

91 CONTINUE 
END 

BLOCK DRTR XYZ 

COMMQN/XNY/NflVG, IFRRM, I BUF 1 ( 1 5 ) , I BUF2 C 5 1 4 ) , 

-BUFSU £ 5 1 4 ) , KOUNT, I REWND, I CNWD, I CNWD1 , 1 CNWD2 , 1 CNWD3 , NPTS, X C 5 1 4 ) 
CQMMON/GRRPH/IGCB£ 192 ) , IPRRM<5 ) , IOBUF (20) , I DCB< 1 44 ) , IFILE (3 ) , 
-XMIN, XMRX, YMIN, YMRX, NCRV,ICRV, VXMIN, VXMRX, VYMIN, VYMRX, NPLT, 

-NFLOT, CRLIB, I£NTRY,MON, I DRY, I YERR, IHOUR, MINUT, IRM 
COMMON ’ T NSTN2 , 1 NSTN3 , 1 NSTN4 , 1 NSTN5 
END 

SUBROUTINE SMOOT £ Y1 , N, IP) 

DIMENSION Y1 £514 ) 

IMRX-N-1 
YI *Y1 £ 1 ) 

DO 5 I-2/IMRX 

IF£IP.EQ.2)YNEW-£YI+2.#Y1 £I)+Y1 £1+1 ))/4. 

IF£IP.EQ.2)YNEW-££11 ./32. >#YI+<1 1 ,/32. )*Y1 £1 )+£ i ,/4 . )*Y1 £1+1 ) 

-+£1 ./IS. )#Y1 £1+2)) 

YI-Y1 £X ) 

YI £1 ) «YNEW 
5 CONTINUE 

IF £IP, EQ. 2 ) YI £N)-£YI+3.#Y1 £N))/4. 

IF£IP.EQ.2)Y1 £N)“££1 ./IB. )»Y1 £1-2 )+£ 1 . /4 . )#Y1 £1*1 ) + (11 ,/32. )*Y1 £1) 
-+£ 1 1 , /32 . )#Y1 £N) ) 

RETURN 

END 

END! 
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tel 




0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 
0010 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
0027. 
0020 
0029 

6655 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0082 
0063 
004S4 


FTN4 , L ORIGINAL PAGE IS 

PROGRAM FNTEMC5) OP POOR QUALITY 

COMMON I YES, NO 

CQMMON/I TI CS/I TI Cl , I TX C2 , 1 TI C3 , 1 TI C4 , 1 TI CS , I TI C6 
COMMQN/XNY/I BUF 1 C 1 5 ) , X C 5 1 2 ) , Y C 5 1 29 , SUM C 5 1 2 ) , NRVG 
COMMON IREWND, ICNWD, ICNWD1 , ICNWD2, ICNWD3 
COMMON I NSTN2 , I NSTN3 , I NSTN4 , I NSTNS 
COMMON XJUNKC10) 

DIMENSION I DAY C 2 ) , I MON C 2 ) , I YEAR C 2 ) , I HOUR ( 2 ) , I MI N C 2 ) 

C 

1 00 FORMRTC //, 8X, " # SHOTS" , 4X, " DATE" , 6X, " TIME" , 6X, " H/V" , 7X, 

T/P" , 7X, " QP/BL" , /, 66 C " »" ) , / ) 

101 F0RMATC10X,Z4,3X,I2,fl1,I2,R1,I2,3X,I2,A1 ,12, 1 X, A2, 4X, Hi , BX, 13, 
-9X,R2) 

103 FORMRTC" SKIPPING tilt ") 

C 

200 FORMRTC/," HOW MANY FILES SHOULD I SKIP?") 

202 FORMRTC" SHOULD I REWIND THE TAPE? C Y/N)?" ) 

68 FORMRT(AI) 

204 FORMRTC/," HOW MANY FILES SHOULD I BACKSPACE?") 

42 FORMRTC" YOU MUST ANSWER YC YES) OR NCNO) ONLY II, TRY AGAIN") 
900 FORMRTC" FILE READ TERMINATED ### TRASH FILE ENCOUNTERED") 

C 

14 WRITE Cl, 202) 

READC 1 , 60 ) IRNS 

IF Cl ANS. EQ. I YES)GQ TO 351 

IF CIANS,EQ.NO)GQ TO 91 

IF C CIRNS. NE . I YES) . AND. CIRNS. NE . NO ) ) WRITE C 1 , 42 ) 

GO TO 1 4 

351 CRLL^EXECCS, IREWND) 

91 WRITE Cl, 200) 

RE AD Cl,*) NSKIP 
IFCNSKIP.EQ.O) GO TO 11 
WRITE Cl ,103) 

WRITEC 1,100) 

NFILE-O 

DO 10 1*1, NSKIP 

CALL EXEC Cl ,ICNWD,IBUF1 ,15) 

96 IF CIBUF1 (8 ) . EQ. 1 )IBUF1 C8)-ITIC3 
IF (I BUF 1 (8 ) . EQ. 2 ) I3UF1 <0)-ITIC4 
IFCIBUF1 (lO).EQ.I )IBUF1 C10)-ITIC5 
IF (IBUF1 CIO). EQ. 2 )IBUF1 C 1 0 ) -I TICS 

WRITEC 1, 101 )IBUF1 C7),IBUF1 Cl ),ITIC2,IBUF1 C2 ) ,ITIC2, IBUF1 (3), 

-I BUF 1 C 4 ) , ITI Cl , I BUF 1 C 5 ) , I BUF 1 ( 6 ) , I BUF 1 ( 8 ) , I BUF 1 C 9 ) , I BUF 1(10) 

I F C NFI LE . NE . 0 ) GO TO 9? 

IFCIBUF1 (D.EQ.13 ) GO TO 99 
CRLL EXECC3, ICNWD1 ) • 

10 CONTINUE 

1 1 CONTINUE 

WRI TEC1 ,204) 

RE AD Cl,*) NBRCK 
IF CNBRCK. EQ. 0) GO TO 13 
NBACK1 « r NBACK+1 
DO 12 1*1, NBRCK1 
CALL EXEC (3, I CNWD2 ) 

12 CONTINUE 

CRLL EXECC3, I CNWD1 ) t 

13 CONTINUE 

DO 14 1-1,512 
XCD-0.0 

14 CONTINUE 
WRITE Cl, 207) 

207 FORMRTC" INPUT THE TIME FOR THE BLOCKED CASE IN THE FORM",/, 


008S 



» XX<HR.>,YY<MIN.> I.E. 10,15") 

0066 



READC1,#)IH0UR<1 ),IMIN<1 ) 

0067 



XJUNK< 1 ) -IHQURC 1 ) 

0066 



XJUNKC2 ) -IMINC 1 ) 

0063 



WRITE <1,2 08) ORIGINAL PAGE 

0070 



OF POOR QUALITY 

0071 




0072 




0073 




0074 


208 

FORMflTC" INPUT THE TIME FOR THE OPEN CASE IN THE FORM",/, 

0075 



" XX<HR. > , YY<MIN, I. E. 10,47") 

0076 



XJUNKC3 ) -IH0URC2 ) 

0077 



XJUNK<4 ) -IMINC2 ) 

0070 



NFILE-2 

0079 



WRITE< 1 ,100) 

0080 



DO 30 K-1 , NFILE 

0081 



CPU- EXEC< 1 , 1 CNWD, IBUF1 ,15) 

0082 



GO TO 96 

0083 


97 

IFCIBUFI < 1 ) .EQ. 1 3 .AND. IBUF1 <4 ) .EQ. 1 3 )G0 TO 99 

0004 



IFCIBUFI <4 ) „ EQ. IH0URC2 ) . AND. IBUF1 <5 ) ,EQ. IMIN<2 ) )G0 TO 16 

0005 



XFCIBUF1C4).EQ.IH0UR<1). AND. IBUF1C5).EQ. IMINC' ))G0 TO 16 

0086 



GO TO 24 

0087 


16 

CONTINUE 

0080 



NAVG-IBUF1 (?) 

0089 



DO 15 1-1,512 

0090 



SUMCI ) -0. 0 

0091 


15 

CONTINUE 

0092 



DO 20 J - 1 , NAVG 

0093 



CALL EXEC < 1 , 1 CNWD, I BUF2 ,514) 

0094 



DO 10 1*3,514 

0095 



SUMCI -2 i "SUMCI -2 )+IBUF2 <1 ) 

0036 


IS 

CONTINUE 

0097 


20 

CONTINUE 

0098 



IFCIBUFI (4). NE.IH0URC2). OR. IBUF1 <5 ) . NE . IMINC2 ) )G0 TO 22 

0099 



IDAYC2 ) -IBUF1 <1 ) 

0100 



IMQNC2) -IBUF1 <2) 

0101 



I YEARC2) -IBUF1 <3 ) 

0102 



CALL VSI4Y C1./NAVG, SUM, 1,Y, 1,512) 

0103 


22 

CONTINUE 

0104 



IFCIBUFI (4). NE.IH0URC1 ). OR. fBUFI C5) . NE. IMINC 1 )>G0 TO 14 

0105 



IDAYC1 )-IBUF1 Cl ) 

0106 



IMONC 1 ) -IBUF1 C2 ) 

0107 



IYEARC 1 ) -IBUF1 C3) 

0108 



CALL VSMY C 1 „ /NAVG, SUM, 1 , X, 1 , 5 1 2 ) 

0109 


24 

CONTINUE 

01 10 



CALL EXECC3,ICNWD1 ) 

01 1 1 


30 

CONTINUE 

0112 



WRITE Cl ,209) 

0113 

C 

209 

. 

0114 

C 



0115 

C 



0116 

C 



0117 

c 



0110 

c 



0119 

c 



0120 

c 



0121 

c 



0122 

c 



0123 

c 



0124 

c 


t 

0125 

c 



0126 

c 



0127 



GO TO 993 

0120 


99 

CONTINUE 

0129 



WRITE <1,900) 

0130 



CALL EXEC C 3,1 CNWD2 ) 27 


27 


0131 

0132 

0133 

0134 

0135 

0136 


CRLL EXECC3, ICNWD1 ) 
STOP 

999 CONTINUE 

GO TO INSTN4 

END 

ENDt 
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0001 

0002 

0003 

0004 

0005 

0006 
000 ? 
0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
001 ? 
0010 

0019 
0C20 
0021 
0022 

0023 

0024 

0025 

0026 
002 ? 

0020 
002 i 

0030 

0031 

0032 

0033 

0034 

0035 

0036 
003? 
0030 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 
004? 
0040 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
0050 

0059 

0060 
0061 
0062 

0063 

0064 


FTN,L 


PROGRAM PTTEM(5 ) 

C'OMMQN/XNY.'NAVG, IFRRM, I BUF 1 ( 1 5 ) # 

-BUFSU ( 5 1 4 ) , KOUNT, I REWMD, I CNWD, I CNWD1 , I CNWD2 , I CNWD3 , NPTS, X ( 5 1 4 ) 
COMMON/GRAPIVI GCB ( 1 92 ) , I PRRM ( 5 ) , I OBUF ( 20 ) , I DCB ( 1 44 ) , I FI LE C 3 ) 

- , XMI N , XMRX , YMI N, YMRX , NCRV , I CRV , VXMI N , VXM AX , VYMI N , VYMRX , NPLT 
-, NPLQT, CALI B, I ENTRY, MON, I DAY, I YEAR, I HOUR , MI NUT, I AM 
COMMON I NSTN2 , I NSTN3 , I NSTN4 , I NSTN5 
DIMENSION DIST(6) 


FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FQRMATC 
FORMAT ( 
FORMATC 


SPO" ) 
SP1 " ) 
SP2 H ) 
SP3" ) 
SP4" ) 
SP5" ) 
SP6" ) 
SP?" ) 
SP0’ ) 
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IENTRY-IENTRY+1 
DO 10 1-1, NPTS 
X(I)-I 
) CONTINUE 
LU“IPRRM( 1 ) 

ID-IPRRMC2) 

IFCID.EQ. 1 . AMD. I ENTRY, EQ. 1 ) 

-CALL PL0TR(IGCB,ID,4,LU,I0BUF,2O) 

IP<ID.EQ. 2. AND. IENTRY.EQ. 1 ) 

-CALL PLOTR C I GCB . ID. 4 , LU, IOBUF ,20) 
IP-MOD(ICRV,B> 

IF (IP.LE . 0) IP-0 

IF (IP .EQ, 1) WRITE (LU, 101) 

IF (IP .EQ. 2) WRITE(LU, 102) 

IF (IP .EQ. 3) WRITE (LU, 103) 

IF (IP .EQ. 4) WRITE (LU, 104) 

IF (IP , EQ. 5) WRITE (LU, 1 05 ) 

IF (IP .EQ. 6) WRITE (LU, 106) 

IF (IP .EO. 7) WRITE (LU, 107) 

IF (IP .EQ. 0) WRITE (LU, 108) 

CALL SETARdGCB, 1.5) 

CALL VI E WP ( I GCB , VXMI N , VXM AX , VYMI N , VYM AX ) 

CALL WI NDW C I GCB , XMI N , XM AX , YMI N , YMAX ) 

HGT-3.0 

AR-0.4 

SANG-0.0 

ICOQR-O 

CALL CSI 2E ( I GCB , HGT , AR , SANG , I CO OR ) 

CALL FXDdGCB, 0) 

XTIC" (XMIN-XMAX)/’S. 

YTIC" ( YMRX -YMI N) ✓’5 . 

IF (IFRRM. NE . 1 ) GO TO 42 
CALL FRAME (I GCB) 

CALL FRAME (I GCB) 

CALL LAXESdGCB, XTIC, YTIC, XMIN, YMIN, 1 . , 1 . ,2. ) 
VXMI 3 - VXMI N- 3. 

VXMR3 ■ VXMRX+6 . 

VYMI 3 -VYMI N- 3. 

VYMR3 » VYMAX+3 . 

DO 43 IDIST-1,6 

DISTdBIST) “(XMIN-dDIST-1 ) #XTIC)*CRLI B*. 3/2 . 
i CONTINUE 

XCENT • ( VXMAX+VXMI N ) /Z . 

YCENT- ( VYMAX+VYMI N ) /Z . 

VXDIF -VXMIN-8 . 

VYDIF -VYMIN-6 . 


0065 

0066 
006 ? 
0066 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 
007 ? 
0070 
0079 
0090 
0081 
0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 
0035 
0096 
009 ? 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 
010 ? 
0100 

0109 

0110 
0111 
0112 

0113 

0114 

0115 

0116 
011 ? 
0110 

0119 

0120 
0121 


CAUL VIEWPCIGCB, VXDIF, VXMAX, VYDIF, VYMRX) 

CALL hll NDW CIGCB, VXDI F , VXMAX , VYDIF, VYMRX ) 

CALL LORGCIGCB, 1 ) 

CALL MOVE ( I GCB,XCENT, VYDIF) 

CALL CPLQTCIGCB, -13.0,0. ,0) 

CALL LABEL Cl GCB) nOtAtm*. 

WRITE CLU, 431 ) C DI ST C I DI ST ) , I DI ST " 1 , 6 ) UnlGINAl PAGE IS 

431 FORMAT C1X,5CF4,1,6X),F4.1 ) °F POOR QUALITY 

C 

VXMA2-VXMA3-3. 

CALL MOVE ( I GCB , VXMA2 , VYMI N ) 

CALL LDIRCIGCB,1 .57) 

IF (IBUF1 (8) .EQ.2) GO TO 30 
CALL LABEL (I GCB) 

WRI TE CLU, 432 ) NAVG, MON, I DRY, I YEAR, I HOUR, MI NUT, I AM 
GO TO 31 

30 CONTINUE 

CALL LABEL (X GCB) 

WRI TE C LU, 433 ) NAVG, MON, I DAY, I YEAR, I HOUR, MI NUT, I AM 

31 CONTINUE 

CALL LDIRCIGCB, 0. ) 

432 FORMAT < 1 4," SHOT AVERAGE - HORIZONTAL" 

-6X, 1 2 , " , 1 2 , " , 1 2 , 4X, 1 2 , " : " , 1 2 , 1 X, R2 ) 

433 FORMAT (I 4," SHOT AVERAGE - VERTICAL 
-6X, 12, " , 12, , 12, 4X, 12, " : " , 12, 1 X, R2 ) 

CALL MO VE ( I GCB , VXMI N , VYMI N ) 

CALL VIEWPCIGCB, VXMIN, VXMAX, VYMIN, VYMRX) 

CALL MI NDW C TGCB , XMI N , XMAX , YMIN, YM AX ) 

C 

C 

42 CONTINUE 

DO 20 J-3,NPTS 
IFCJ.GT.3) GO TO 21 

I F < I CR V . NE . 1 ) CALL MOVECIGCB,XCJ),BUFSUCJ) ) 

IFCICRV.EQ. 1 )CRLL MOVE CIGCB, XCJ ) , BUFSUC J ) ) 

GO TO 20 

21 I F C I CR V . NE . 1 ) CALL DRRWCIGCB, XCJ ), BUFSUC J ) ) 

IFCICRV.EQ. 1 )CRLL DRRWCIGCB,X( J),BUFSUCJ) ) 

20 CONTINUE 

XTI ME ■ XMI N+3 . 2 # C XMAX - XMI N)^5 . 

YTI ME • YMAX-I CRV* C YMRX-YMI N ) ^1 0 . 

CALL MOVE Cl GCB, XTI ME, YTI ME) 

CALL LABEL Cl GCB) 

IFCICRV.EQ. 1 ) WRITE CLU, 201 ) 

C IFCICRV.EQ. 2) WRITE CLU, 202 ) 

IFCICRV.EQ. 2) WRITE CLU, 203 ) 

201 FQRMRTC" AVE . NO FIT") 

C 202 FORMAT C" AVE. 3-PTS. FIT") 

203 FORMAT C" AVE. 5-PTS. FIT") 

CALL XMI TCI GCB ) 

41 CONTINUE 

I F C NPLT . EQ . NFLO T . AND . I CRV . E Q . NCRV ) WRI TE C LU , 1 00 ) 

IF CNPLT. EQ. NPLOT. AND. ICRV. EQ.NCRV) CALL PLOTRCIGCB, ID, 0) 

IFCICRV.EQ. 1 ) GO TO INSTN3 

IF CICRV. NE . 1 ) GO TO INSTN2 

END 

ENDS 


30 
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FTN L OMGINAL PAGE 

* PROGRAM DPLOT OF POOR QUALITY 

COMMON/XNY/X (512),Y(512),NP 

COMM&N/GRAPH/I GCB (192), I PRfiM <5 ) , IQBUF (20) , I DCB (144),IFILE(3) 

• , XMI N , XMAX, YMIN, YMRX, NCRV, ICRV, VXMIN, VXMAX, VYMIN, VYMRX, NPLT 
NFLOT 

COMMON/CRED/X JUNK ( 20 ) 

COMMON INSTN1 , INSTN2, INSTN3, INSTN4 
DIMENSION VPNQM(4, 4 ) 

I NTEGER REDUC ( 3 ) , PTRED (3 ) , DAFITO) 

EQUIVRLENCE (IPRRM(1 ),LU), C I PRAM ( 2 ) , I D ) 

DRYA VPNOM/1 0 . , 70. , 57 . , 97 . , 1 0. , 70. , 7 . , 47 . , 85 . , 1 45 . , 57 . , 97 . , 

-85. ,145. ,7. ,47./ 

DRTR REDUC/2HRE , 2HDU, 2HC / 

DATA PTRED/2HPT , 2HRE , 2HD / 

DATA DAFI T/2HDA, 2HFI , 2H2 / 

C 

C THIS CALL PICKS UP THE LU AND ID OF THE PLOTTING DEVICE 
CALL RMPAR(IPRRM) 

C 

600 FORMAT (" THIS PRORGAM IS THE MAIN PROGRAM FOR SETTING UP PLOTS",/, 

" IT READS IN THE PLOT DEVICE LU AND ID THRU RMPRR. .„.",/, 

" I.E. RU,GPLOT, 1,1 PLOTS ON THE CRT OR",/, 

" RU, GPLOT ,32,2 PLOTS ON THE PLOTTER",/) 

601 FORMAT (" QUESTIONS WILL BE ASKED CONCERNING HOW MANY PLOTS",/, 

" SCALING OF EACH (INDIVIDUALLY) RND PEN NUMBER",/) 

WRITE (1,600) 

WRITE (1,601 ) 

C 

C THIS ASSURES THAT A VALID LU/ID COMBINATION IS USED 
IF(LU.L£.0 .OR. ID.LE.O) GO TO 90 
IF(ID.EQ.I) LU-1 

IF(ID.EQ.2. AND. (LU.NE.32. AND.LU.NE.33) ) GO TO 90 
IF(ID. GT. 2) GO TO 90 


C 

C 


C 


C 


c 


300 FORMAT( 1 X, " HOW MANY FRAMES DO YOU WANT?") 

400 FORMAT (" INPUT VIEWPORT COORD.: VXMIN, VXMAX, VYMIN, VYMAX" , /, 

FOR FRAME", 13," DEFAULT VALUES: ", 3 (F3 . 0, ","), F3 . 0) 

401 FORMAT (" INPUT WINDOW COORD. XMI N, XMAX, YMIN, YMRX" , /, 

DEFAULT VALUES: 0,10,0,1000") 

402 FORMAT (" HOW MANY CURVES ON FRAME" , 12, " ?" ) 


502 FORMAT (//" 

_ ii 


ii 


***#****#**************************************'• 
* PLACE PEN IN HOLDER *" 

***********************************************" 


ASSIGN 101 TO INSTN1' 

ASSIGN 102 TO INSTN2 
ASSIGN 103 TO INSTN3 
ASSIGN 104 TO INSTN4 

WRITE (1,300) 

READd,*) NPLOT 
DO 92 NPLT- 1, NPLOT 
VXMIN -0. 

VXMRX-O. 

KNOM-MQD(NPLT, 4 ) 

IF (KNQM . EQ. 0) KNOM-4 

WRITE ( 1 , 400) NPLT, ( VPNGM( JNOM, KNOM) , JNOM-1 , 4 ) 
READd,#) VXMIN, VXMAX, VYMIN, VYMRX 
IF (VXMIN. NE . 0. .OR. VXMAX. NE.O.) GO TO 30 
VXMIN- VPNOMd ,KNOM) 

VXMAX -VPNQM(2, KNOM) 


' > 
//) 


31 
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0066 
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0090 
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0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
01 1 1 
0112 
0113 
01 1*4 

0115 

0116 
0117 


VYMIN»VPN0M(3,KN0M) 

VYMRX"VPNQM(4, KNQM) 

30 CONTINUE 
C 

XMIN-O. 

XMAX-O. 

WRITE (1,401 ) 

RERD< 1 , #) XMI N , XMRX, YMI N, YMAX 
IFCXMIN.NE.O. .OR. XMRX.NE.O.) GO TO 31 
XMIN-O. 

XMRX-10. 

YMIN-Q. 

YMRX-1000. 

31 CONTINUE 


ORIGINAL PAGE EG 
OF POOR QUALITY 


C 

C WRITE (1,402) NPLT 

C RERD( 1 , *) NCRV 

C DO 92 ICRV-1,NCRV 

NCRV" 2 
ICRV-1 

IF (LU. EQ. 32 ) WRITE ( 1 , 502 ) 
C 


CRLL EXEC(B, REDUC) 
101 CONTINUE 


NP-XMRX/( . 15*XJUNK(5) ) 
IF (NP. GT. 500) NP-500 
CALL EXEC(8, PTRED) 

102 CONTINUE 


CRLL EXEC(a.nRFIT) 
103 CONTINUE 


ICRV-2 

CRLL EXEC <8, PTRED) 
104 CONTINUE 


92 CONTINUE 
GO TO 91 

90 CONTINUE 
WRITE (1,999) 

999 FQRMRT( IX, "IMPROPER LU OR ID... USE PARAMETERS IN RU COMMAND") 

91 CONTINUE 
END 

BLOCK DRTR XY2 
COMMON/'XNY^X ( 5 1 2 ) , Y(5 1 2 ) , NP 

COMMON/GR RPH/I GCB (1923,1 PRAM ( 5 ) , 1 0 BUF ( 20 ) , I DCB (144),IFILE(3> 
-, XMIN, XMRX, YMI N, YMRX, NCRV, ICRV, VXMIN, VXMAX, VYMIN, VYMAX, NPLT 
-,NPU3T • 

COMMON^CRED/X JUNK ( 20 ) 

COMMON I NSTN1 , I NSTN2 , I NSTN3 , I NSTN4 

END 

END* 
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PROGRAM REDUCC5 > 0F POOR QUALITY 

COMMON /X NY/ XC51 2 ) , YC51 2 ) , NP 
COMMON /CRED/ XJUNKC20) 

COMMON INSTN1 , INSTN2,INSTN3,INSTN4 
DIMENSION I BUFFI C 1 5 > , IBUFF2 C51 4 ) , SUMC5 1 2 ) 

DIMENSION I DAY <2 ) , I MON <2 ) , I YE ARC 2 ) , I HOUR C 2 ) , IMINC2 ) 


C 


100 FORMAT C20X, " SHOTS H/V T/P OP/BL" ) 

101 F0RMRTC2X,X2,"/",I2,"/",I2,2X,I2," :",X2, IX, A2,3I5,3X, A2,5I5) 

102 FQRMATC2X,2I5) 

103 FORMAT C» SKIPPING " ) 

C 

200 FORMAT (" HOW MANY FILES SHOULD I SKIP?") 

201 FORMATC" HOW MANY FILES HAVE TO BE READ TO OBTAIN THOSE FILES?") 

202 FORMATC" SHOULD I REWIND THE TAPE? CY OR N)") 

203. FORMAT CR1 ) 

204 FORMATC" HOW MANY FILES SHOULD I BACKSPACE?" ) 

205 FORMATC" INPUT THE LU FOR THE OUTPUT") 

206 FORMATC" IT IS ASSUMED THAT THE SKIP AND BACKSPACE COMMANDS" ,/, 

" POSITIONED THE TAPE CLOSE TO THE DESIRED DATA. NOW",/) 

207 FORMATC/" INPUT THE TIME FOR THE BLOCKED CASE IN THE FORM",/, 

~ " 9 45" ) 

200 FORMATC/" INPUT THE TIME FOR THE OPEN CASE IN THE FORM",/, 

" 9,46") 

209 FORMATC/" INPUT THE TYPE OF SMOOTHING THAT YOU WANT DONE",/, 

" 0 - NO SMOOTHING",/," .1 - RECTANGULAR SMOOTHING",/, 

" 2 - TRIANGULRR SMOOTHING CDIST WILL BE FULL WIDTH AT” 

- ' ," HALF MAX)",/," 3 - GAUSSIAN CDIST IS 2X 1-SIGMA WIDTH)") 

210 FORMATC/" INPUT 'DIST', THE SIZE IN KILOMETERS OF THE SMOOTHING" 

," FILTER.") 

C 

ICNWD-10B 

IREWND-410B 

ICBK-1410B 

ICFWD-1310B 

IYES-1HY 

N0-1HN 

C 

WRITE Cl, 202) 

RE AD Cl, 203) IANS 

IF CIANS.EQ.IYES) CALL EXECC3,IREWND) 

C 

WRITE Cl ,200) 

READC1,*) NSKIP 
IF CNSKIP. EQ. 0) GO TO 11 
WRITE Cl, 103) 

WRITE Cl, 100) 

DO 10 1-1 , NSKIP 

CALL EXECC 1 , ICNWD, IBUFF1 ,15) 

WRITE C 1 , 1 01 ) CIBUFF1 CJ), J-1 ,15) 

IF Cl BUFFI C 1 ) . EQ. 13 ) GO TO 99 
CALL EXECC3,ICFWD) 

10 CONTINUE 

11 CONTINUE 
C 

WRITE Cl, 204) 

RE AD Cl,*) NBACK 

IF CNBACK. EQ, 0) GO TO 1 3 i 

NBACK1 -NBACK+1 
DO 12 I «1 , NBACK 1 
CALL EXEC C 3,1 CBK) 

12 CONTINUE . „ 

CALL EXECC3, ICFWD) ^ 


0065 


13 

CONTINUE 



0066 

C 





006? 



DO 14 1-1, 5T3 



0068 



X(I i-0.0 



0069 


14 

CONTINUE 



0070 

C 





0071 



WRITE (1,206) 



0C*»2 



WRITE (1 , 207) 



0073 



READ( 1 , #) I HOUR ( 1 ) , IMIN( 1 ) 



0074 



XJUNK (1 )-IH0UR(1 ) ■- ^ 



0075 



XJUNK(2) -IMIN( 1 ) ORIGINAL PAGE fS? 



0076 

c 


OF POOR QUALITY 



0077 



WRITE (1,200) 



0078 



RERD( 1 , #) I HOUR (2), I NI N(2) 



0079 



XJUNKO ) -IHQUR(2 ) 



0080 



X JUNK (4)-lMlN(2) 



0081 

c 





0082 



WRITE (1,201 ) 



0083 



READd,#) NFILE 



0084 

c 





0005 



WRITE (1,100) 



0006 



DO 30 K-1, NFILE 



008? 



CRLL EXEC (1,1 CNWD, I BUFF 1,15) 



0080 



WRITE (1,101) (I BUFFI (I), 1-1, 15) 



0089 



IF (I BUFFI (1 ) . £Q. 13 .RND. I BUFFI (4 ) . EQ. 1 3 ) GO TO 99 



0090 



IF (I BUFFI (4) . EQ.IH0UR(2) .RND. IBUFF1 (5 ) . EQ. IMINC2 ) ) GO 

TO 

16 

0091 



IF (I BUFFI (4) . EQ.IHOURd ) . RND. TBUFF1 (5 ) . EQ. IMINd ) ) GO 

TO 

16 

0092 



GO TO 24 



0093 


16 

CONTINUE 



0094 



NRVG-IBUFF1 (7) 



0095 



XJUNKd 7 ) -NRVG 



0096 



DO 15 1-1,512 



0097 



SUM (I ) “0. 0 



0090 


15 

CONTINUE 



0099 

c 





0100 



DO 20 J-1 , NRVG 



0101 



CRLL EXEC(1 , XCNWD, IBUFF2,514) 



0102 



DO 10 1-3,514 



0103 



SUM (I -2 ) -SUM(I -2 )+IBUFF2 (I ) 



0104 


10 

CONTINUE 



0105 


20 

CONTINUE 



0106 

c 





0107 



YRVG-1 .0/FL0AT(NRVG) 



0108 



IF (I BUFFI (4) .NE. I HOUR (2) .OR. IBUFF1 (5) .NE. IMIN(2)) GO 

TO 

22 

0109 



IDRY(2)-IBUFF1 (1 ) 



0110 



IM0N(2)-IBUFF1 (2) 



0111 



I YERR(2 ) -I BUFF 1 (3) 



0112 



XJUNK ( 6 ) -I DRY ( 2 > 



0113 



XJUNK ( 7 ) -I MON ( 2 ) 



0114 



XJUNK ( 8 ) - I YEAR ( 2 ) 



0115 



XJUNK ( 9 ) -I BUFF 1 (6) 



0116 



XJUNKC 16) -IBUFF1 (8)' 



0117 



CRLIB-FLORTd BUFFI (9) )/10. 



0118 



CRLL VSMYCYAVG,SUM, 1 , Y, 1 ,512) 



0119 


22 

CONTINUE 



0120 



IF (I BUFF 1 (4 ) . NE . IH0URC1 ) .OR. IBUFF1 (5 ) . NE . IMINd ) ) GO TO 24 

0121 



I DRY ( 1 ) -I BUFF 1(1) 



0122 



IMONd ) “IBUFF1 (2) 



0123 



I YEAR (1 ) -I BUFFI (3) 



0124 



CRLL VSMY(YRVG,SUM,1,X, 1,512) 



0125 


24 

CONTINUE 



0126 

c 





012? 



CRLL EX£C(3, XCFWD) 



0120 

c 


: 



0129 


30 

CONTINUE 



0130 



WRITE (1,209) 




VMM 


0131 


RERD(1,#) IREC 


0132 


XJUNK( 1 0 ) -IREC 


0133 


WRITEd ,210) 


0134 


RERD(1,#) DIST 


0135 


IFCIREC.NE.O) CRLL SliOTHCY, DIST, GRUB, IREC) 


0136 


IF (IREC. NE . 0) CRLL SMOTH(X, DIST, CRLIB, IREC ) 


0137 


XJUMK( 1 9 ) -DIST 


0139 

0139 


CRLL VSUB (Y,1,X,1,Y, 1,512) 
XJUNK(5 ) -CRLIB 

np'S mL 
OF POOR QUALITY 

0140 


DO 40 1*1,512 

0141 


X(I ) -I#CRLIB*. 1 5 

0142 

0143 

C 

40 CONTINUE 


0144 


GO TO 999 


0145 


99 CONTINUE 


0146 


WRITEd , 900) 


0147 


CRLL EX£C(3, ICBK) 


0140 


CRLL EXECC3, ICFWD) 


0149 


900 FQRMRTC FILE RERD TERMINATED ### TRRSH FILE 

ENCOUNTERED" ) 

0150 


999 CONTINUE 


0151 


GO TO INSTN1 


0152 


END 


0153 


SUBROUTINE SMOTH(Y, DIST, CRLIB, IREC) 


0154 


DIMENSION Y(512),T(512),W(801 >,WSUMC512) 


0155 

C 

IREC-1 RECTRNGULRR SMOOTHING 


0156 

C 

IREC-2 TRI RNGULRR SMOOTHING 


0157 

C 

IREC-3 GRUSSIRN SMOOTHING 


0150 


NP-DTST/(CRLIB#.15) 


0159 


IF (IREC. NE . 1 ) NP»2.0#DIST/(CRLIB#. 15) 


0160 


IFCIREC.EG.3) NSIG-O. 5#DI ST/ (CRLIB#. 1 5 ) 


Oiei 


IF (IREC. EU. 3) MU-Z#NSIG 


0162 


IF (IREC. EQ. 3 ) NP-2#MU+1 


0163 


IF (NP. GT. 800) NP-000 


0164 


NP2-NP/2 


0165 


NP-NP2#2+1 


0166 


DI ST "NP#CRLI B# . 1 5 


0167 


IF (IREC. NE . 1 ) DI ST -NP#CRLIB#. 15/2.0 


0160 

C 

IFCIREC.EQ.3) DIST-NP#CRLIB». 1 5/2. 0 


0168 

0170 

c 

DO 10 1-1, NP2 


0171 


IF (IREC. EQ. 1 ) W(I)-1. 


0172 


IF (IREC. EQ. 2 ) W(I )»FLORT(I )/FL0RT(NP2+1 ) 


0173 


IF (IREC. EQ. 3 ) W(I ) -EXP( - . 5# (FLOAT (I -MU- 1 )/FLORT(NSIG) )**2) 

0174 


W(NP+1-I)-WCI) 


0175 


10 CONTINUE 


0176 


W(NP2+1 )-1 .0 


0177 


DO 15 1-1,512 


0178 


MIN-I-NP2 


0179 


MRX-I+NP2 


0100 


IF(MIN.LT.1 )’ MIN-1 


0181 


IF (MRX. GT. 51 2 ) MRX-512 


0182 


T(I ) -0, 0 


0183 


WSUMd ) -0. 0 


0104 


DO 15 J -MIN, MRX 


0185 


INDEX-J-I+NP2+1 


0196 


T(I ) »T(I )+W( J-I+NP2+1 ) #Y ( J ) 


0107 


WSUMd) -WSUM(I )+W( J-I+NP2+1 ) 


0180 


15 CONTINUE 


0109 


DO 20 1-1,512 


0190 


Yd)-T(I) /WSUM d ) 


0191 


20 CONTINUE 


0192 


RETURN 


0193 


END 


0194 


END* 
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PROGRAM DRFI2 (5 ) OF POOR QUALITY 

COMMON /XNY/ XC512 ) , SCSI 2 ) ,NP 
COMMON /CRED/ XJUNKC20) 

COMMON I NSTN 1 , I NSTN2 , I N5TN3 , I NSTN4 

DIMENSION PARA<4),PSTRC3,4),SIGC3),SXGTC4),TEMPC4>,IC<4),JIN<4) 

C 

200 FORMAT <" INPUT THE MIN DIST. AND MAX DIST. FOR THE FIT, EG 0,6") 

201 FORMATC* INPUT THE ORDER IN WHICH YOU WANT TO FIT. EG. 2, 1,3, 4",. 

" WILL FIT ARER, KAPPA, ALPHA, FOCUS.") 

CALIB-XJUNKC5) 

WRITE (1 ,200) 

READ ( 1 , *) RMIN,RMRX 
N1 "RMIN/C . 1 5*CRLIB) 

N2 “RMRX/ C . 1 5*CRLI B ) 

WRITE<1 ,201 > 

READ< 1 , *) JINC1 >,JIN<2),JINC3), JINC4) 

C 

C 

RCEN«(RMAX+RMIN)/2. 

11- RCEN/(.15*CALIB) 

12- RMIN/(.15*CRLIB) 

X3-RMRX/C.15*CALIB> 

G1 ■< 1 . O-RCEN/2 . )**2 
G2»<1 .0-RMIN^2. )**2 
G3»<1 .0-RMAX/2. >**2 

T1 -EXP<-.3*RCEN)/SGRTCCRCEN*1000. )**2+< .040*1 . 0E+05/1 . 06 )**2*G1 ) 
T2-EXP<“.3*RMIN)/SQRTC CRMIN*1000. )**2+< .040*1 .0E+05/1 .06)**2*G2> 
T3"EXPf-.3*RMAX)/SQRT( (RMAX*1000. )**2+C .040*1 .0E+05/1 .06)**2*G3) 
PSTR< 1 , 1 > -SCI 1 )/T1 
PSTR<2,1 )-S(I2)/T2 
PSTR<3, 1 ) "S£I3 )/'T3 
C 

PSTR< 1,2) -0.04 
PSTR(2,2)-0.025 
PSTR<3,2)-0.055 
C 

PSTRC1 ,3) -0.3 
PSTR<2,3)-0. 1 
PSTR(3,3)-0.5 
C 

PSTR< 1 , 4 ) -2. 

PSTR<2, 4 ) ■ .5 
PSTR(3, 4 ) » 5. 

C 

DO 93 N-1, 10 
DO 13 J-1,4 
1F( J.EQ* 1 ) KJ-JINC1 ) 

IFCJ.EQ.2) KJ-JINC2) 

IF ( J . EQ. 3 ) KJ-JINC3) 

IF <J.EQ.4) KJ-JINC4) 

C 

C KJ « 1 KAPPA 

C 2 AREA 

C 3 ALPHA 

C 4 FOCRL RANGE 

C 

DO 12 1-1,3 
C 

DO 11 L-1 , 4 
PRRA(L) »PSTR( 1 , L) 

11 CONTINUE 

PARA < K J ) -PSTR ( I , KJ ) 

C 
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0094 
QQ95 

ooii 

0097 

0098 
0039 
0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 
01 1 1 
0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 
0128 

0127 

0128 

0129 

0130 


CALL SI GMR ( PRRR ( 1 ) , PARR ( 2 ) , PRRR ( 3 > , PRRR ( 4 ) , SI G C I) , S , N 1 , N2 , X ) 
12 CONTINUE 

T1-(SIG(1 )-SIG(2))/(PSTR(1 , KJ ) -PSTR(2, KJ ) ) 
T2-(SIG(2)-SXG(3))/(PSTR(2,KJ)-PSTR(3,KJ>) 
R2-(T1-T2)/'(PSTR(1,KJ)-PSTR(3,KJ)) 
R1-T2-R2*(PSTR(2,KJ)+PSTR(3,KJ)) 

RO-SI G C 3 ) -fll #PSTR ( 3 , KJ ) -R2 #PSTR ( 3 , KJ ) **2 
UPDRT--R1/(2.*R2) 


14 

15 


20 


21 


I MAX • 1 

I F ( SI G ( 2 ) . GT , SI G C 1 ) ) IMRX-2 
IF(SIG(3).GT.SIG(IMAX)) IMRX-3 
I MIN* 1 

IF(SIG(2).LT.SIG(1 )) IMIN-2 
IF(SIG(3).LT.SIG(IMIN)> IMIN-3 
ICEN-1 

IF CICEN.EQ. IMRX .OR. ICEN. EQ. IMIN) ICEN-2 
IF (ICEN.EQ.IMRX .OR. ICEN. EQ. IMIN) ICEN-3 
SLOPE-(PSTR(IMIN,KJ)-PSTR(ICEN,KJ>)/(SIG(IMIN)-SIGdCEN)) 

PRRR (KJ ) "PSTR (I MI N, KJ ) -SLOPE#SI G (I MI N ) 

I F C PRRR ( K J ) . LT . 0 . ) PRRR < K J ) -RBS ( PRRR (KJ ) ) 

IFLRG-0 

I F ( UPDRT . LT . PSTR (IMRX, KJ ) . AND . UPDRT . GT . PSTR (I MI N , K J ) . AND . R2 . GT . 0 . ) 
IFLRG-1 

I F ( UPDRT . LT . PSTR ( I MRX, K J ) . AND . UPDRT . GT. PSTR (I MIN, KJ ) . RND. R2 . GT. 0. ) 

PRRR(KJ) •UPDRT 
IF2-1 

WRITE (31 , 1Q1 ) PflRfl(KJ),PSTR(IMIN,KJ>,SIG(lMlN),PSTR(ICEN,KJ), 
SIG(ICEN) , PSTR (IMRX, KJ ) , SXG (IMRX) 

F0RMRT(E1 2 . 3 , 6E3 .3 ) 

IF (PRRR(KJ ) . GT. 1 0. #PSTR(IMIN, KJ ) ) PRRR ( KJ ) -PSTR (IMIN, KJ ) 

DO 14 1-1,10 

IF (I . GT. 1 ) PRRfl (KJ ) ■ (PSTR ( IMIN, KJ )+PRRR (KJ ) )/ / 2 . 

CALL SIGMR(PRRR(1 ) ,PRRfl(2 ) , PRRRO ) ,PRRR<4 , SDEV, S, N1 ,N2,X) 

IF (SDEV.LTi SIG(IMIN) ) GO TO 15 

CONTINUE 

IF2-0 

CONTINUE 

TEMP ( 1 ) "PSTR ( 1 , KJ ) 

TEMP ( 2 ) “PSTR ( 2 , KJ ) 

TEMP ( 3 ) "PSTR ( 3 , KJ ) 

TEMP ( 4 ) «PRRR ( KJ ) 

SIGT( 1 ) ■SIG( 1 ) 

SIGT(2)-SIG(2) 

SIGTC3) -SIGC3) 

SIGT(4 ) "SDEV 
YSTOR-SIGTd ) 

ICC 1 ) -1 
DO 21 1-1,4 
DO 20 L»I,4 

I F ( SI GT ( L ) . GT . YSTQR ) GO TO 20 
YSTOR-SIGTd) 

IC(I )-L 
CONTINUE 
INDX-IC(I ) 

SI GT ( I NDX ) -SI GT ( I ) 

SIGT(I ) -YSTQR 
PSTOR-TEMPdNDX) 

TEMP (I NDX) -TEMP (I) 

TEMPd ) -PSTOR 
CONTINUE 


ORIGINAL PAGE m 
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PSTR (1,KJ)- TEMPd ) 
PSTR(2,KJ)-TEMP(2) 
PSTR ( 3 , K J ) °TEMP ( 3 ) 


37 


i 0131 

013?. 

c 

16 

CONTINUE 


0133 



WRITE (1,1 00) SIGTC1 >, CPSTR< 1,1) ,1*1,4), IFLAG, 

IF2 

, 0134 

0133 

c 

100 

FORMAT C" SDEV, XK, A, RL, F -» ,569.3,212) 


0136 


13 

CONTINUE 


013? 

0138 

0139 

0140 

c 

99 

CONTINUE 

XJUNK (10) -PSTR (1,3) 
XJUNKM1 )-PSTR(1 ,2) 

ORIGINAL PAGE 
OF POOR QUALITY 

0141 



XJUNK (12) "PSTR C 1 , 4 ) 


0142 



XJUNK (13) *PSTR (1,1) 


0143 



XJUNK ( 1 4 ) -RMIN 


'• 01 44 

014b 

c 


XJUNK < 1 5 ) ■RMAX 

. 

■ 0146 



CALL CALC ( PSTR < 1 , 1 ) , PSTRC 1,2), PSTRC 1,3), PSTR< 1 , 4 ) , X, S ) 

; oh? 



GO TO INSTN3 


0148 
1 0149 

c 


END 


0150 

0151 

c 


SUBROUTI NE SI GMA ( XK, A, AL, F , SDE V, S, N1 , N2 , X ) 


0152 



DIMENSION S(512),X(512) 


0153 



3UM-0. 


0154 



SUM2-0. 


0155 



DO 10 I-N1,N2 


; 0156 



T1-XK*EXPC-RL#Xa>> 


, , 0157 



T2-C1 .0-X(I)/F)*#2 


y r 0158 



T2-T2#1 ,0E+10*A*R/<1 .06*1 .06) 


y is 0159 



T2-<X(I)#1000. )*#2+Y2 


l. L CA 

» - - V 1 W V 



Cr*DI (*«T< yOADT t TO * 

WWI IkNV — 1 « *■ W%RM | I 1M r 


■ :• oi6i 



SUM-SUM+ C S (I ) -SCALC ) #*2 


i 0162 



SUM2 -SUM2+ < S (I ) -SCRLC ) 


* ' 0163 


10 

CONTINUE 


0164 



SDEV«SQRT< SUM/FLOAT CN2-N1+1 )) 


0165 



SUM2-SUM2/FL0ATCN2-N1+1 ) 


0166 

c 


WRI TE < 3 1 , 1 00 ) SDE V , SUM2 , XK, R , AL, F 


0167 

0166 

c 

100 

FORMATS SDEV,SUM2,XK,A,AL,F«" ,6E9.3) 


0,169 

0170 



RETURN . 

END 


0171 

c 




0172 

0173 

c 


SUBROUTINE CALCCXK, A, AL, F, X, Y) 


0174 



DIMENSION X(512) , YC512) 


0175 



DO 10 1-1,512 


0176 



T1 -XKXEXPC -AL#X(I ) ) 


0177 



T2- ( 1 .0-XCI)^F>*#2 


[ 0170 



T2-T2*A#A*1 .0E+10/C1 .06*1 .06) 


i 01 79 



T2-T2+CXCI >*1000. >**2 


0180 



Y<I)»TVSQRTCT2> 


5 0181 


10 

CONTINUE 


0182 



RETURN 


i 0183 



END 


0104 



END* 


■ i 



t 



&PTRED T-00004 IS ON CRQ0Q12 USING 00010 BLKS R-0000 


0001 

0002 

0003 

0004 

0005 

0006 
000 ? 
0000 

0009 

0010 
0011 
0012 

0013 

0014 
00'* 5 
0016 
001? 
0010 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 
0035 
003? 

0038 

0039 

0040 

0041 

0042 

0043 
004 4 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
0050 
0059 
OOSO 
0061 
0062 

0063 

0064 


FTN, L 


100 

101 

102 

103 

104 

105 

106 
107 
100 
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c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM PTRED(5 > 

COMMON/XNY/X.< 512), ¥(512), NP 

COMMON/GRRPH/I GCB < 1 92 ) , I PRRM C 5 ) , 1 0 BUF ( 20 ) , I DCB ( 1 44 ) , I FI LE ( 3 ) 
*, XMIN, XMAX, YMXN, YMRX, NCRV, ICRV, VXMIN, VXMRX, VYMIN, VYMRX, NPLT 
•,NPLQT 

COMMON/'CRED/'XJ UNK ( 20 ) 

COMMON INSTN1 , INSTN2, INSTN3, INSTN4 
FORMAT (" SPO" ) 

FORMRTC'SPI" ) 

FORMAT <”SP2" ) 

FORMAT (" SP3" ) 

FORMAT <"SP4" ) 

FORMAT CSP5" ) 

FORMAT <" SP6" ) 

FORMAT <» SP7" ) 

FORMATC" SP0" > 

LU-IPRRMC1 ) 

ID-IPRAMC2) 

IF (ID. EQ. 1 .AMD. NPLT. EQ. 1 .AND. ICRV. EQ. 1 ) 

■CALL PLOTRf IGCB, ID, 4, LU, IQBUF, 20> 

IF(ID.EQ.2.AND.NPlT.EQ.1 . AND. ICRV. EQ.1 ) 

•CALL PLQTRdGCB, ID, 4, LU, IQBUF, 20) 

IP«MOD(ICRV,0) 

IF (IP. LE . 0) IP-8 
IF(IP .EQ. 1) WRITE (LU, TQ1 ) 

2) WRITE (LU, 102) 

3) WRITE (LU, 103) 

4) WRITECLU, 104 ) 

5) WRITE (LU, 105) 

6) WRITE (LU, 106) 

7) WRITE (LU, 10?) 

0) WRITE(LU, 100) 

1 . AND.LU.NE.32) WRITE (LU, 101 ) 


IF(IP 

IF(IP 

IF (IP 

IF dP 

IFdP 

IFdP 

IFdP 

IFdP 

IFdP 

IFdP 

IFdP 

IFdP 

IFdP 

IFdP 

IFdP 

CALL 

CALL 


. E Q . 
.EQ. 
.EQ. 
.EQ. 
.EQ. 
.EQ. 
.EQ. 
.EQ. 
.EQ. 

H LQ. 
.EQ. 
. EQ. 
. EQ • 
.EQ. 
. EQ. 


2. AND.LU.NE.32) 

3. AND.LU.NE.32) 

4. AND.LU.NE.32) 

5. AND.LU.NE.32) 

6 . AND. LU. NE . 32 ) 

7. AND.LU.NE.32) 


WRITE (LU, 102) 
WRITE (LU, 103 ) 
WRITE (LU, 104) 
WRITE (LU, 105) 
WRITE (LU, 106) 
WRITE (LU, 107) 
WRITE (LU, 1 08 ) 


8. AND.LU.NE.32) 

SE^ARdGCB, 1.5) 

VI E WP ( I GCB, VXMI N , VXMRX , VYMIN, VYMflX ) 


CRLL WINDWd GCB , XMIN, XMHX , YMI N , YMRX ) 

HGT-3.0 

AR-O. 4 


SRNG-0.0 

ICOOR-O 

CRLL CSI 2E ( I GCB, HGT, RR, SANG, I COOR ) 

CRLL FXDdGCB, 0) 

XTIC«(XMIN-XMAX)/S. 

YTI C-( YMRX- YMI N)/5. 

IFdP. EQ.1) CRLL FRAME (IGCB ) 

IFdP. EQ.1) CALL FRRME (IGCB) 

IF (IP . EQ . 1 )CALL LRXE S ( I GCB , XTI C , YTI C , XMI N , YMIN, 1 . , 1 . , 2 . ) 

DO 20 J - 1 , NP 

IF( J . GT. 1 ) GO TO 21 

CALL MOVE (IGCB, X(J),Y(J)) 

GO TO 20 

21 CALL DRRW(IGCB,X(J),Y( J) ) 

20 CONTINUE 

IF (ICRV. GT. 1 ) GO TO 10* 

VXMI3-VXMIN-3. 

VXMR3 ■ VXMRX+6 . 


39 


0065 

0066 
0067 
0060 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0001 
0082 

0003 

0004 

0005 
0086 
0087 

. 0080 
1 0009 

; 0090 

: oo9i 

1 0092 

? 0093 

| 0094 

' 0095 

i 0096 
0097 
0090 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 
0107 
0100 

0109 

0110 
01 1 1 
0112 

0113 

0114 

0115 
' 0116 
j 011? 

! 0110 
* 0119 

: 0120 

0121 
0122 
0' 73 
i C. A 
j 0125 
? 0126 
1 01 27 

I 0128 

ft 0129 
c 0130 


VYMI3-VYMIN-3. 

VYMA3-VYMRX+3. 

CRLL VI E WP (I GCB , VXMI 3 , VXMR3 , VYMI 3 , V YM A3 ) 
CRLL WI NDW ( I GCB , VXMI 3 , VXMR3 , VYMI 3 , VYMR3 ) 
CRLL L0RG(IGCB,1) 

VXMR2-VXMR3-3. 

BIST «XJUNK( 1 9 ) 

NRVG-XJUNKC17) 

IM“XJUNK(6 ) OF 

IDRY®XJUNK(7> • OF 

I Y£R**XJUNK( 8 ) 

IREC®XJUNK(10) 

ILOOK-XJIJNK(IS) 

CRLL MOVE ( I GCB , VXMR2., VYMI N ) 

CRLL LDXRCIGCB, 1.57) 

IF (XLOQK. EQ. 2 )G0 TO 30 
CRLL LABEL (IGCB) 

WRITE (LU, 432 )NRVG, IM, IDRY, IYER 

432 FORMAT ( 1 4, " SHOTS-RVERRGE HORIZHNTRL" , /, 
-12X,I2," /",I2," /" ,12) 

GO TO 31 

30 CONTINUE 

CRLL LRBELdGCB) 

WRITE (LU, 433 )NRVG,IM, IDRY, IYER 

433 FORMAT (1 4," SHOTS- RVERRGE VERTICAL",/, 
-1 2X, 12, " /" ,12,"/" ,12) 

31 CONTINUE 

CRLL LBIRdGCB, 0. ) 

XCENT ■ ( VXMRX+ VXMI N)/2 . 

YCENT® ( VYMflX+VYMI N ) '2 . 

VXDIF-VXMIN-2.#HGT 

vydif-vymih-2.*hgt 

CALL VI E WP < I GCB, VXDI F , VXMRX , VYDI F , VYMRX ) 
CRLL WI NDW ( I GCB , VXDI F , VXMRX , VYDI F , VYMRX ) 
CRLL LQRGdGCB, 1 ) 

CRLL MOVE (IGCB, XCENT, VYDIF) 

CRLL CPLOT ( I GCB ,-13.,Q. ,0) 

CRLL LABEL <1 GCB) 

WRITE (LU, 700) 

CRLL MOVE (IGCB, VXDIF, YCENT) 

CRLL CPLOT (I GCB, 0. , -5. ,0) 

CRLL LDIRCIGCB, 1.57) 

CALL LABEL (I GCB) 

WRITE (LU, 701 ) 

CRLL LDIR(IGCB, 0. ) 

• CRLL VIEWPdGCB, VXMI N, VXMRX, VYMI N, VYMRX) 
CRLL WI NDW ( I GCB , XMI N , XMRX> YMI N , YMRX ) 

XTI ME « XMI N+3 . 2 * < XMflX - XMI N ) /5 . 

YTIME "YMRX- (YMRX- YMI N) /I 0. 

CALL MOVE (IGCB, XTI ME, YTIME) 

IH»XJUNK(1 ) 

IM-XJUNI<(2) 

CALL LABEL (I GCB) 

WRITE (LU, 702) IH,IM 
YTIME-YMRX-2.*CYMRX-YMIN)/10. 

CRLL MOVE (IGCB, XTI ME, YTIME) 

IH»XJUNK(3) 

IM’“XJUNK(4 ) 

CRLL LABEL d GCB) 

WRITE (LU, 703 ) IH,IM 

XTI ME ® XMI N+2 . 6 * ( XMRX-XMI N ) /5 . 

YTIME® YMRX - 3 . * ( YMRX - YMI N ) / 1 0 . 


K * 

kvi 

BsfS 

f-jl 
£. | 
?• *1 


0131 CALL MOVE ( I GCB, XTI ME , YTI ME ) 

0132 GO TO (505,506,507,509), (IREC+1 ) 

0133 SOS CONTINUE 

0134 CRLL LABEL (I GCB) 

0135 WRITE (LU, 705) 

0136 705 FORMATC NO -SMOOTHING" ) 

0137 GO TO 9 

0138 506 CONTINUE 

0139 CRLL LRBELCIGCB) 

0140 WRITE (LU, 706) 

0141 706 FORMAT (" RECTANGULAR-SMOOTHING" ) 

0142 GO TO 9 

0143 507 CONTINUE 

0144 CRLL LABEL(IGCS) 

0145 WRITE(LU,707) 

0146 707 FORMAT (" TRIANGULAR-SMOOTHING" ) 

0147 GO TO 9 

0140 500 CONTINUE 

0149 CRLL LABEL(IGCB) 

0150 WRITE (LU, 708) 

0151 708 FORMAT (" GRUSSI AN-SMOOTHING" ) 

0152 9 YTIME-YMRX-4 . *( YMRX-YMIN)/1 0. 

0153 XTI ME - XMI N+2 . 7 * ( XM AX - XMI N ) /5 . 

0154 CRLL MOVE ( I GCB, XTI ME , YTI ME ) 

0155 CALL LABEL (I GCB) 

0156 WRITE (LU, 709) BIST 

0157 709 FORMATC FILTER WIDTH-" , F4 . 2, " KM.") 

0158 10 CONTINUE 

0159 IFdCRV.EQ. 1 )• GO TO 25 

0160 RMRX-XJUNK(IS) 

0161 RMIN-XJUNK( 1 4 ) 

0162 AL-XJUNK(IO) 

0163 R-XJUNK(II) 

0164 F -XJUNK( 1 2 ) 

0165 XK-XJUNK( 1 3 ) 

0166 XTIME-XMIN+3 . 2*(XMAX-XMIN)/5 . 0 

0167 YTIME»YMRX-5.*(YMAX-YMIN)/10. 

0168 CRLL MOVE ( I GCB, XTI ME , YTI ME ) 

0169 CALL LABEL (I GCB) 

0170 WRITE CLU, 704) RL, A,F,XK 

0171 XG01 -XMIN+RMIN 

0172 XG02-XMIN+RMRX 

0173 YG01 -YMAX-(YMRX-YMIN)/10. 

0174 CRLL MOVEdGCB, XG01 ,YMRX) 

0175 CALL PLOTdGCB, XG01 , YMAX, -2 ) 

0176 CRLL PLOTdGCB, XG01 ,YG01 , -1 ) 

0177 CALL MOVE (I GCB, XG02, YMRX) 

0178 CRLL PLOTdGCB, XG02, YMAX, -2) 

0179 CALL PLOTdGCB, XG02,YG01 , -1 ) 

0190 25 CONTINUE 

0181 700 FQRMRTC" DISTANCE, KILOMETERS ") 

0182 701 FORMATC’ RELATIVE SIGNAL, ARB. UNITS") 

0183 702 FORMAT ("BLOCKED" ,13," ,12) 

0104 703 FORMRTC' OPEN 

0185 704 FORMATC' AL-" , El 0. 3, /, " A •" , El 0. 3, /, " F 

0186 CALL XMITCIGCB) 

0187 IF (NPLT, EQ. NPLOT. AND. ICRV. EQ. NCRV) 

0188 IF (NPLT. EQ. NPLOT. AND. ICRV. EQ. NCRV) 

0109 IFdCRV.EQ. 1 ) GO TO INSTN2 

0190 IFdCRV.EQ. 2) GO TO INSTN4 

0191 END 

0192 END* 


ORIGINAL PAGE IS 
OF POOR QUALITY 


/ 


,£10.3,/," XK-" , E10. 3) 


WRITE (LU, 100) 

CRLL PLOTRCIGCB, ID, 0) 


41 


